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Preface 

The  study  involved  in  this  thesis  lias  been  an  attempt  to  find 
out  whether  there  was  some  possible  application  of  the  Mossbauer 
Effect  to  accelerometers  used  in  inertial  navigation  systems.  After 
one  year  of  research  in  this  area  I  can  be  neither  comprehensive  nor 
conclusive  in  the  presentation  of  this  problem.  In  writing  the  thesis 
it  has  been  my  hope  that  the  material  presented  will  serve  the  follow¬ 
ing  purposes: 

First,  to  provide  the  applications  engineer  with  a  succinct, 
but  not  oversimplified,  introduction  to  the  Mossbauer  Effect. 

Second,  to  outline  clearly  some  of  the  major  limitations  in 
the  application  of  the  Mossbauer  Effect  and  the  fundamental  reasons 
for  these  limitations. 

Third,  to  bring  into  focus  by  means  of  an  order  of  magni¬ 
tude  analysis  some  of  the  more  promising  areas  of  investigation  which 
have  come  to  my  attention. 

At  one  point  in  the  study  I  considered  introducing  classified 
performance  figures  from  the  more  recent  advanced  techniques  in 
inertial  navigation.  This  would  have  provided  a  more  precise  basis 
for  comparing  the  suggested  techniques  to  the  present  state  of  the 
art.  However,  in  view  of  the  generality  of  my  presentation,  I  con¬ 
cluded  that  the  inclusion  of  the  classified  information  would  not  add 
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so  much  to  the  Thesis  as  to  justify  the  concomitant  administrative 
burden. 

The  area  of  investigation  was  suggested  by  Mr.  Paul  Polishuk 
of  the  Flight  Control  Laboratory,  Aeronautical  Systems  Division.  Mr. 
Polishuk  was  most  helpful  in  providing  reference  material  and  in 
periodically  reviewing  the  study.  Dr.  Bernard  Kaplan  acted  as  my 
thesis  adviser,  providing  periodic  counseling  and  review  of  my  efforts. 
Finally,  I  am  grateful  to  my  wife,  for  typing  my  notes  as  well 

as  assembling  and  typing  the  final  manuscript. 


•Jack  E.  Hesse,  Jr. 
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Abstract 

An  analytical  investigation  of  the  application  of  the  Mossbauer 
Effect  to  accelerometers  is  presented.  The  Mossbauer  Effect  is 
the  emission  and  absorption  of  gamma  photons  without  energy  losses 
due  to  the  recoil  of  the  emitting  and  absorbing  nuclei.  Resonant 
absorption  so  obtained  is  sensitive  to  changes  in  photon  energy  of 
one  part  in  10^  -  10^.  Fundamental  limitations  in  the  application 
of  this  phenomenon  are  the  need  for  statistical  counting  of  photons 
and  the  uncertainty  in  the  energy  of  the  nuclear  gamma  decay.  An 
order  of  magnitude  analysis  indicates  that  the  isotope  Fe^,  which 
has  favorable  characteristics  such  as  a  relatively  high  cross  section 
for  resonant  absorption  and  a  relatively  high  ratio  of  photon  energy 
to  uncertainty  in  photon  energy  (3  x  1012),  is  the  best  available 
isotope.  The  Mossbauer  Effect  is  most  sensitive  to  a  shift  in  energy 
due  to  relative  velocity  between  the  source  and  absorber,  allowing 
sensing  of  velocities  of  10"^  -  10"^  cm/sec.  Techniques  are  presented 
and  analyzed  that  use  the  latter  phenomenon  to  measure  angular  and 
linear  velocity  and  angular  position  of  instrument  components.  A 
typical  estimate  for  one  second  of  counting  shows  that  angular  velo¬ 
city  may  be  measured  to  within  10~G  radians/ sec  and  angular  position 
to  within  2  seconds  of  arc.  An  application  of  these  techniques  to 
selected  accelerometers  does  not  indicate  a  promising  increase  in 
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accuracy.  Uncertainties  in  parameters  not  related  to  the  Mossbauer 
Effect  such  as  friction,  hysterisis,  mass  unbalance,  and  angular 
measurement  limit  the  possible  accuracy  to  about  one  part  in  10^. 

A  promising  application  is  found  in  the  Rotating  Pendulum  Accelero¬ 
meter  which  offers  the  possibility  of  a  low  threshold  sensitivity,  10"^  g 
or  less. 
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AN  INVESTIGATION  OF  THE  APPLICATION 
OF  THE  MOSSBAUER  EFFECT 
TO  ACCELEROMETERS  USED  FOR  INERTIAL  NAVIGATION 

I.  Introduction 

This  thesis  presents  the  results  of  an  investigation  to  deter¬ 
mine  whether  there  are  possible  applications  of  the  Mossbauer  Effect 
to  accelerometers  used  in  inertial  navigation  systems.  The  study  was 
motivated  by  the  following  factors: 

First,  the  approaching  era  of  space  travel  will  place  re¬ 
quirements  of  extreme  accuracy  on  navigation  systems.  This  is 
particularly  true  for  inertial  components  such  as  accelerometers. 

The  Mossbauer  Effect,  although  discovered  in  1958,  has 
already  found  wide  application  in  scientific  experiments  requiring 
extreme  sensitivity  and  precision.  It  is  then  natural  to  question 
whether  the  Mossbauer  Effect  can  be  used  to  advantage  in  other 
applications  requiring  unusual  accuracy,  such  as  the  measurement 
of  acceleration  during  space  missions. 

In  this  chapter  the  above  two  statements  will  be  expanded  in 
some  detail,  and  the  plan  of  the  thesis  will  be  outlined. 

Requirements  of  Accelerometers  Used  for  Ine rtial  Navigation 

The  requirements  for  extreme  accuracy  in  space  navigation 
stem  from  the  great  distances  and  length  of  time  involved  in  the 


1 


GA/Phys/62-7 

flights.  A  relatively  small  error  expressed  as  a  fraction  of  the 
total  distance  from,  say  Earth  to  Mars,  will  represent  an  error 
of  many  thousands  of  miles.  There  are  several  approaches  to  the 
problem  of  space  navigation,  including  tracking  and  controlling  the 
vehicle  from  an  earth  based  command  station,  celestial  navigation, 
and  inertial  navigation.  Future  space  flights  may  employ  any  one  of 
these  methods,  or  more  likely,  a  combination  of  all  three.  Tracking 
and  controlling  from  an  earth  station  saves  weight  in  the  assembly 
of  navigation  components  carried  by  the  vehicle;  however,  when  one 
considers  that  it  will  take  several  minutes  for  an  electronic  signal 
to  reach  earth  from  a  vehicle  in  the  vicinity  of  one  of  the  nearer 
planets,  the  disadvantages  of  this  method  become  apparent.  Using 
celestial  navigation,  a  system  will  require  several  "fixes"  taken  over 
a  considerable  interval  of  distance  and  time  in  order  to  determine 
important  flight  parameters  such  as  velocity.  Thus  it  will  be  diffi¬ 
cult  to  monitor  and  control  the  characteristically  short  bursts  of 
powered  flight  by  the  use  of  celestial  navigation.  Therefore,  it  may 
be  necessary,  or  at  least  desirable,  to  revert  to  an  inertial  system 
for  some  part  of  the  total  space  navigation  problem.  The  chief  limi¬ 
tation  in  the  use  of  inertial  systems  for  space  flight  is  the  requirement 
for  extreme  accuracy. 

An  inertial  system  determines  the  position  of  a  vehicle  in 
space  by  twice  integrating  the  accelerations  incident  upon  the  vehicle. 
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The  measurement  of  these  accelerations  is  accomplished  by  the  use 
of  accelerometers.  The  requirement  for  extreme  accuracy  is  illus¬ 
trated  by  the  following: 

Let 

S  =  JJAdt  (1) 

where  S  is  the  position  of  the  vehicle,  A  is  the  acceleration  of  the 
vehicle,  and  dt  is  a  differential  element  of  time.  If  aA  is  a  con¬ 
stant  error  present  in  the  measurement  of  acceleration,  then  AS 
will  be  the  resultant  error  in  position,  given  by 

AS  =  A  A  t2 

T  (2) 

The  significance  of  the  t2  term  for  space  flights  of  long  duration  can 
be  readily  appreciated. 

In  addition  to  a  requirement  for  extreme  accuracy,  it  is  desir¬ 
able  for  purposes  of  space  navigation  to  have  accelerometers  that 
measure  very  small  values  of  acceleration.  One  reason  for  this 
requirement  is  the  advent  of  small  thrust  propulsion  devices  such  as 
ion  engines  and  solar  sails.  For  example,  consider  a  device  which 
produces  10"2  lbs  of  thrust  on  a  1000  lb  vehicle.  The  resultant 
acceleration  will  be  10-D  times  the  value  of  acceleration  due  to  gravity 
at  the  earth's  surface.  Another  source  of  motivation  for  developing 
sensitive  accelerometers  is  the  possible  use  of  gravity  gradient  tech¬ 
niques  for  space  navigation.  These  methods  ensue  from  the  fact  that 
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even  in  "free  fall"  only  the  center  of  mass  of  a  vehicle  is  subject 
to  an  exact  weightless  condition.  Positions  of  the  vehicle  at  a  greater 
or  smaller  distance  from  the  attracting  body  than  the  distance  from 
that  body  to  the  center  of  mass  of  the  vehicle  are  subject  to  small 
values  of  specific  force.  These  small  specific  forces  are  due  to  the 
gradient  of  gravity,  which  acts  in  accordance  with  the  inverse  square 
law.  A  description  of  some  possible  gravity  gradient  techniques  is 
contained  in  Appendix  A. 

Applications  of  the  Mossbauer  Effect 

A  complete  description  of  the  Mossbauer  Effect  will  be  given 
in  Chapter  II.  In  this  section  it  will  suffice  to  say  that  the  Mossbauer 
Effect  is  a  form  of  nuclear  resonance  fluorescence  of  gamma  rays. 

It  is  distinguished  from  other  forms  of  nuclear  resonance  fluorescence 
in  that  gamma  rays  are  emitted  and  absorbed  without  recoil  of  the 
individual  emitting  and  absorbing  nuclei. 

Rudolf  L.  Mossbauer,  a  young  German  Physicist,  discovered 
the  effect  named  after  him  in  1958  (Ref  16).  His  discovery  rapidly 
proved  to  be  an  invaluable  tool  for  scientific  research.  Using  the 
Mossbauer  Effect,  R.  V.  Pound  and  G.A.  Rebka,  Jr.  were  able  to 
measure  the  effect  of  gravity  on  the  energy  of  gamma  photons  "fal¬ 
ling"  through  laboratory  distances.  This  is  a  direct  confirmation  of 
Einstein’s  predicted  red  shift  of  spectral  lines  in  a  gravitational  field 
(Ref  19,20).  The  Mossbauer  Effect  has  been  used  to  study  lattice 
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vibrations  (Ref  30),  nuclear  magnetic  effects  (Ref  4,6),  and  lifetimes 
of  excited  nuclei  (Ref  15:20).  As  will  be  demonstrated  in  Chapter  II, 
by  using  the  Mossbauer  Effect  one  can  detect  relative  velocities  of 

n  c 

the  orders  10 to  10“3  cm/sec.  For  comparison,  the  end  of  the 
minute  hand  on  a  watch  moves  at  about  2x10“^  cm/sec. 

In  view  of  the  above,  one  would  not  be  surprised  to  see  the 
Mossbauer  Effect  go  into  space  along  with  the  MASER*  and  the  LASER, 
which,  like  the  Mossbauer  Effect,  exhibit  sharply  resonant  phenomena. 
However,  in  contrast  to  the  LASER  and  the  MASER,  little  or  no  work 
has  been  done  to  realize  an  engineering  application  from  the  Mossbauer 
Effect.  Therefore,  this  study  must  start  with  the  theory  and  scienti¬ 
fic  applications  of  the  Mossbauer  Effect  and  proceed  in  a  very  general 
and  tortuous  way  to  seek  an  engineering  application  in  the  field  of 
accelerometers  used  for  inertial  navigation. 

Thesis  Plan 

Chapter  n  provides  a  review  of  the  theory  and  major  experi¬ 
ments  involved  in  the  Mossbauer  Effect.  The  most  important  isotopes 
are  discussed,  and  relationships  showing  the  effect  of  different  physical 
phenomena  on  the  Mossbauer  Effect  are  developed.  In  Chapter  III  the 
thesis  problem  is  initiated  by  outlining  and  explaining  general  limitations 

*  MASER  or  LASER,  respectively  Microwave  or  Light  Amplification 
by  Stimulated  Emission  of  Radiation. 
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applicable  to  any  experiment  performed  with  the  Mossbauer  Effect. 
Chapter  IV  contains  a  description  and  an  order  of  magnitude  analysis 
of  several  hypothetical  schemes  for  using  the  Mossbauer  Effect  to 
measure  angular  position  and  angular  or  linear  velocity  of  an  instru¬ 
ment  component.  In  Chapter  V  the  various  types  of  accelerometers 
are  classified,  and  selected  accelerometers  are  analyzed  with  respect 
to  an  application  of  the  Mossbauer  Effect.  Chapter  VI  contains  the 
conclusions  and  recommendations  derived  from  the  study. 
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II.  A  Brief  Review  of  the  Mossbauer  Effect 

As  stated  previously,  the  Mossbauer  Effect  is  a  form  of  nuc¬ 
lear  resonance  fluorescence  of  gamma  rays.  Nuclear  resonance 
fluorescence  implies  the  excitation  of  a  particular  energy  transition 
in  one  nucleus  by  a  source  of  energy  resulting  from  an  identical 
energy  transition  in  another  nucleus.  In  the  Mossbauer  Effect, 
nuclear  resonance  fluorescence  is  accomplished  by  emission, 
absorption,  and  reemission  of  gamma  rays.  Further,  these  gamma 
rays  are  emitted  and  absorbed  without  recoil  of  the  emitting  and 
absorbing  nuclei. 

This  chapter  contains  the  explanation  and  significance  of  the 
Mossbauer  Effect.  The  basic  theory  is  presented;  then  a  brief  re¬ 
view  of  the  major  experiments  involved  in  the  development  of  the 
Mossbauer  Effect.  Two  of  the  most  significant  Mossbauer  isotopes 
are  discussed.  Finally,  relationships  showing  the  effect  of  differ¬ 
ent  physical  phenomena  on  the  Mossbauer  Effect  are  developed. 
Theory 

The  theory  of  the  Mossbauer  Effect  consists  of  an  understand¬ 
ing  of  first,  why  energy  losses  due  to  recoil  of  a  nucleus  occur 
during  the  emission  or  absorption  of  a  gamma  ray,  and  second,  how 
these  recoil  losses  can  be  eliminated. 

Consider  a  free  nucleus,  at  rest.  If  the  nucleus  undergoes 
an  energy  transition,  E^,  and  thereby  emits  a  gamma  photon,  the 
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the  law  of  conservation  of  momentum  demands  that  part  of  the 
energy  of  transition  be  given  up  to  the  recoil  of  the  nucleus.  This 
energy  lost  to  recoil  is  denoted  R,  given  by 

R  =  P2  (3) 

where  P  is  the  momentum,  and  M  is  the  mass  of  the  nucleus.  By 

momentum  conservation, 

P  =  P  =  Ey_  (4) 

C 

where  p  is  the  photon  momentum,  is  the  photon  energy,  and  C 
is  photon  velocity.  By  substituting  (4)  into  (3),  the  following  results: 

R=EV2  (5) 

mat 

Since  energy  is  conserved,  Ey  =  Et  -  R.  A  similar  argument  can 
be  made  to  show  that  for  absorption  to  take  place,  the  energy  re¬ 
quired  of  the  photon  is  given  by  E^  =E^  +  R. 

The  value  of  Ej.  is  not  an  exact  quantity,  but  is  indeterminate 
in  the  amount  r  ,  given  by  the  Heisenberg  Uncertainty  Principle 

r  =  h  (6) 

Tm 

where  ft  is  Planck's  constant  divided  by  2n,  and  Tm  is  the  mean 
lifetime  of  the  excited  state  of  the  nucleus  (Ref  29:113).  For  nuclear 
transitions,  R  >>  r,  therefore  the  energy  separation  between  the 
emitted  photon  and  the  energy  required  for  resonant  absorption  is 
approximately  2R,  as  illustrated  in  Figure  1.  If  this  were  the  only 
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consideration,  then  the  energy  difference  2R  could  be  supplied  by 
some  artificial  means  such  as  by  moving  the  source,  thereby  im¬ 
parting  a  Doppler  shift  in  the  energy  of  the  emitted  photons.  How¬ 
ever,  in  nature  the  nuclei  involved  are  not  at  rest  as  implied  by 
the  previous  model,  but  are  moving  with  respect  to  one  another. 

In  a  gas  the  atoms  containing  the  nuclei  are  in  random  thermal 
motion  and  in  a  solid  the  atoms  vibrate  about  points  in  a  crystal 
lattice.  These  motions  induce  a  broadening  in  the  energy  distribu¬ 
tion  of  the  emitted  photons,  rendering  artificial  compensation  for 
recoil  losses  a  difficult  task. 

The  essential  feature  of  Mossbauer's  discovery  was  that, 
under  certain  conditions,  atoms  in  a  solid  emit  or  absorb  gamma 
photons  without  recoil  of  the  individual  nuclei.  It  is  possible  to 
appreciate  how  this  happens  from  fundamental  results  of  quantum 
mechanics.  The  atoms  in  a  solid  oscillate  with  discrete  energies, 
given  by  E  =  h  v  (Ref  10:122).  E  is  the  energy  of  the  oscillators 
and  v  is  the  frequency  of  oscillation.  The  allowed  energies  are 
discrete,  therefore  only  certain  specific  transitions  are  permitted. 
Suppose,  then  that  all  the  atoms  in  a  crystal  are  in  state  Ej  and 
the  next  allowed  state  is  E£.  Suppose  further  that  the  transition 
energy,  E2  -  Ej,  is  twice  the  value  of  R,  the  recoil  energy  from 
a  gamma  photon  emission.  If  a  number  of  atoms,  N,  emit  photons, 
and  N/2  atoms  undergo  the  transition  from  state  1  to  state  2,  then 
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the  total  energy  alloted  from  recoil  of  all  N  atoms  will  be  account¬ 
ed  for,  since  E2  -  Ej  =2R.  This  will  leave  N/2  atoms  to  emit 
photons  without  recoil.  Momentum  is  conserved  because  eventually, 
the  recoil  from  a  single  event  is  absorbed  by  the  entire  crystal. 

Application  of  the  preceeding  to  a  more  realistic  model  of  a 
crystal  makes  the  mathematics  considerably  more  complicated,  but 
the  physical  principles  are  the  same.  The  emission  of  a  photon 
without  recoil  is  referred  to  as  a  zero  phonon  process.  A  phonon 
is  the  quantum  of  energy  (or  frequency)  involved  in  the  transition 
of  an  atom  in  a  lattice  from  one  state  to  another  state.  The 
fraction,  f,  of  zero  phonon  processes  that  occur  is  approximated 
by 

f  =exp  3  R _  (  1+  2n2T2)  ~/  (7) 

2  k  0  d  (  3od2  ) 

where  k  is  Boltzman's  constant,  T  is  the  ambient  temperature  of 
the  crystal,  and  0d  is  the  Debye  Temperature.  This  is  a  good 
approximation  where  T  <<  Od  (Ref  20:440). 

The  parameter  0d  is  defined  by 

°d  -  hv  max  (8) 

k 

where  v  max  is  the  maximum  frequency  at  which  any  atom  oscillates 
in  the  crystal  lattice  (Ref  5:39-45).  With  this  definition  in  mind, 
one  can  appreciate  the  following  requirements  for  obtaining  a  detect¬ 
able  Mossbauer  Effect.  A  high  0d  is  necessary  in  order  that  there 
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be  lattice  transitions  whose  energy  is  larger  than  the  free  recoil 
energy  of  the  nucleus.  The  requirement,  as  seen  from  (7),  is  that 
R  <  <k  Gd*  Another  requirement  is  that  T  <<  G^.  Intuitively,  a 
low  T  will  mean  that  the  atoms  in  the  crystal  are  in  the  lower 
energy  states  where  the  transition  energies  are  larger.  These  re¬ 
quirements  are  diagrammed  in  Figure  2. 

If,  in  fact,  the  above  requirements  are  met,  then  a  signifi¬ 
cant  fraction  of  the  gamma  photons  can  be  emitted  or  absorbed  by 
zero  phonon  processes.  The  cross  section,  cr ,  for  resonant  absorp¬ 
tion  is  given  by 

o  -  o_q  f  f  (9) 

2 

where  f  is  the  fraction  of  zero  phonon  occurrences  in  the  emitting 
lattice,  f'  is  the  fraction  of  such  occurrences  in  the  absorber  lattice, 
and  oQ  is  a  basic  cross  section  for  electromagnetic  absorption, 
given  by 

an  =  \  2  2Ie  +  1  1  (10) 

2  n  2lg  +  1  1  +  q 

X  is  the  gamma  photon  wavelength,  Ie  is  the  nuclear  spin  in  the 

excited  state,  Ig  is  the  nuclear  spin  in  the  ground  state,  and  a  is 

the  internal  conversion  coefficient  (Ref  30:208,  4:222,  8:7).  The 
above  expressions  assume  an  isotopic  abundance  of  100%  of  appro¬ 
priate  isotopes  in  the  absorber  lattice.  If,  for  example,  only  2%  of 
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the  atoms  in  the  absorber  lattice  were  the  proper  isotope  for  Moss- 
bauer  absorption,  then  the  actual  cross  section  per  atom  of  the  en¬ 
tire  absorber  would  be  only  2%  of  the  value  given  by  (9). 

If  the  Mossbauer  a  is  plotted  against  energy  as  shown  in 
Figure  3,  a  Breit-Wigner  curve  of  width  2  r  results.  The  width  is 
just  twice  the  value  given  by  the  Heisenberg  Uncertainty  Principle 
for  the  nuclear  transition  energy  and  results  from  an  overlap  of 
the  two  curves  shown  in  Figure  1  (Ref  14:6).  It  is  precisely  the 
narrow  width  of  this  curve  that  makes  the  Mossbauer  Effect  a  deli¬ 
cate  and  sensitive  tool  for  scientific  experiments.  Numerical 
examples  will  be  discussed  shortly,  but  first  a  short  resume  of  the 
important  experimental  milestones  in  the  history  of  the  Mossbauer 
Effect  is  given. 

Experiments 

Prior  to  Mossbauer' s  discovery  attempts  were  made  (with 
some  success)  to  observe  nuclear  resonance  fluorescence  by  com¬ 
pensating  for  the  recoil  losses  in  various  ways.  Mossbauer  himself 
was  engaged  in  just  such  an  attempt  when  he  discovered  the  phenom¬ 
ena  of  recoilless  emission  and  absorption.  He  was  observing  res¬ 
onance  fluorescence  of  the  129  -  kev  gamma  ray  from  Ir1^.  For 
this  transition,  R  is  0.05  ev,  and  the  broadening  of  the  transition 
energy  due  to  thermal  motion  at  room  temperature  is  about  0. 1  ev. 
At  room  temperature  the  emission  and  absorption  spectrum  overlap 
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considerably,  therefore  resonance  fluorescence  can  be  observed. 

In  order  to  reduce  this  scattering  by  reducing  the  overlap,  Moss- 
bauer  cooled  both  source  and  absorber.  He  expected  that  the 
decrease  in  Doppler  broadening  of  the  transition  energy  would  lead 
to  a  decrease  in  the  resonance  fluorescence.  Instead,  the  very 
opposite  happened!  (Ref  8:11-12) 

Mossbauer  determined  the  nature  of  his  result  and  recognized 
its  significance,  that  by  eliminating  recoil  losses  he  had  removed 
the  major  cause  of  line  broadening.  He  explained  his  observations 
by  adapting  a  theory  by  W. E.  Lamb,  Jr.  Lamb's  theory  was  develop¬ 
ed  to  explain  the  effect  of  lattice  binding  on  the  capture  of  slow 
neutrons  (Ref  11:190).  Mossbauer's  experiments  were  repeated  by 
other  laboratories,  particularly  at  Los  Alamos  and  Argonne.  All 
of  these  experiments  were  performed  with  Ir^l  and  were  compli¬ 
cated  by  low-temperature  requirements  and  by  the  low  fraction  of 
zero  phonon  processes  that  occur  in  Ir*^  (about  5%  at  0°K).  (Ref 
8:12). 

The  Mossbauer  Effect  assumed  major  importance  for  practi¬ 
cal  applications  when  the  phenomena  was  observed  in  Fe^.  This 
discovery  was  made  independently  at  Harvard,  Harwell,  the  Uni¬ 
versity  of  Illinois,  and  Argonne.  Because  of  its  high  (about 
450°K),  Fe5?  exhibits  a  large  fraction  of  zero  phonon  processes 
at  room  temperature  (about  70%).  In  addition,  the  linewidth,  r, 


13 


GA/Phys/62-7 

involved  in  the  Fe^  Mossbauer  Effect  is  very  narrow,  allowing  very 
precise  experiments  (Ref  8:12-13). 

A  typical  Mossbauer  experiment  is  shown  schematically  in 
Figure  4.  Photons  from  a  source  are  passed  through  an  opening 
in  shielding  and  allowed  to  strike  the  absorber,  usually  a  foil  only 
a  few  millimeters  thick.  In  the  case  of  an  absorption  experiment 
the  photons  that  pass  through  the  foil  are  counted  by  a  detector  plac¬ 
ed  directly  behind  the  foil.  If  conditions  such  as  temperature  and 
relative  motion  are  properly  controlled,  Mossbauer  absorption  will 
take  place  and  the  counting  rate  will  be  a  minimum.  If  conditions 
are  changed  such  as  to  reduce  the  probability  of  Mossbauer  absorp¬ 
tion,  the  counting  rate  will  increase.  It  is  true  that  the  absorbed 
photons  will  be  reradiated,  but  the  reradiation  will  take  place  iso¬ 
tropically.  By  proper  arrangement  of  the  geometry  of  the  experi¬ 
ment  this  effect  can  be  minimized.  On  the  other  hand,  if  the  de¬ 
tector  is  placed  as  shown  in  Figure  4  for  a  scattering  experiment, 
only  the  reradiated  photons  will  be  counted. 

In  order  to  fully  realize  the  significance  of  the  quantities 
involved  in  such  an  experiment,  the  properties  of  some  important 
Mossbauer  isotopes  will  be  discussed. 

Isotopes 

A  list  of  isotopes  for  which  the  Mossbauer  Effect  has  been 
observed  is  given  in  Appendix  B.  Two  of  these  isotopes  will  be 
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Schematic  Diagram  of  a  Typical 
Mossbauer  Experiment 
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discussed,  Fe57  which  has  the  best  combination  of  desirable 
Mossbauer  properties,  and  Zn67  which  exhibits  the  smallest 
linewidth,  r. 

The  following  decay  scheme  is  of  interest  for  the  Fe57 
Mossbauer  Effect. 


c  57 
'-'0 

— > 

(Fe57)*  +  p  + 

(ID 

(Fe57)* 

— > 

(Fe57)**  +  v  (123  kev) 

(12) 

(Fe57)** 

— > 

Fe57  +  y  (14. 4  kev) 

(13) 

(Fe57)*  decays  by  (12)  and  (13)  in  91%  of  the  instances.  The  other 

9%  of  the  events  result  in  a  direct  conversion  to  ground  state  by 

emission  of  a  137  kev  gamma  ray.  The  phenomenon  of  interest  is 

(13),  which  happens  without  recoil  in  a  very  large  fraction  of  the  cases. 

In  addition  to  a  high  Qd  and  relatively  low  value  of  R,  Fe57 

57 

has  the  desirable  property  that  its  parent  nucleus,  CQ  ,  decays  with 
a  relatively  long  half  life  of  270  days.  Most  important,  the  mean 
lifetime  of  the  excited  state  in  (13)  is  comparatively  long  for  a  Moss- 
baucr  transition.  Tm  =  1.44  x  10-7  sec.  By  (6),  r=  4.  6  x  10-9  ev. 

With  the  exception  of  Zn®7,  Fe57  exhibits  a  smaller  linewidth  than 
any  other  Mossbauer  isotope.  Pound  and  Rebka,  among  others,  have 
experimentally  observed  this  small  linewidth.  Reasonable  care  in 
the  preparation  of  source  and  absorber  will  yield  a  linewidth  within 

*  or  **  Excited  states  of  the  nucleus. 
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twice  the  value  predicted  by  (6).  (Ref  21:555). 

A  parameter,  Q,  can  be  defined  which  is  a  useful  measure 
of  the  degree  of  resonance  in  a  Mossbauer  experiment. 

Q  =  Ey_  (14) 

r 

For  Fe^7,  Q  =  3  x  10^.  The  interpretation  is  that  any  phenomena 
which  changes  the  energy  of  the  gamma  photon  by  an  order  of  one 
part  in  10^  win  significantly  alter  the  resonant  absorption.  This 
is  seen  from  the  plot  on  a  versus  E  in  Figure  3. 

A  comparison  of  the  properties  of  Zn ®7  with  those  of  Fe®7, 
as  listed  in  Apprendix  B,  will  reveal  that  Fe®7  is  generally  super¬ 
ior  except  for  the  fact  that  Zn®7  has  a  longer  mean  lifetime  of 
the  excited  nucleus,  resulting  in  a  very  small  r.  It  is  difficult  to 
perform  a  Mossbauer  experiment  with  Zn®7.  The  low  0^  of  this 
isotope,  combined  with  a  large  value  of  R,  results  in  a  small  frac¬ 
tion  of  zero  phonon  occurences,  even  at  temperatures  near  0°K. 
However,  the  high  value  of  Q  (about  2  x  10  ^)  that  can  be  obtained 
promises  scientific  applications  of  the  utmost  sensitivity.  A  group 
led  by  P. P.  Craig  has  performed  one  of  the  few  successful  series 
of  experiments  with  Zn®7.  This  group  observed  a  linewidth  about 
twice  the  value  of  r  predicted  by  (6).  (Ref  4:561-564) 

The  usefulness  of  a  small  linewidth  and  high  Q  value  is  made 
apparent  by  the  following  discussion  of  the  effects  of  physical  phe¬ 
nomena  on  Mossbauer  experiments. 
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Physical  Phenomena 


Some  of  the  major  physical  phenomena  which  can  be  observed 
through  the  Mossbauer  Effect  are  relative  velocity,  acceleration, 
temperature  effects,  chemical  shifts,  and  magnetic  effects. 

Relative  Velocity.  The  first  mentioned  of  these  phenomena 
acts  through  the  Doppler  shift,  whereby  the  frequency  of  a  photon 
emitted  by  a  moving  source  is  seen  by  a  stationary  observer  to  be 

v=  vQ  C  •  (15) 

The  source  is  moving  toward  the  observer  with  velocity  v.  v  is  the 
observed  frequency,  and  vQ  is  the  emitted  frequency  (Ref  26:442-444). 
By  subtracting  vQ  from  both  sides  of  the  expression,  the  following 
can  be  obtained: 

Av  =  v  -  Vn=  Vn  V  (16) 

C-v 

If  V«  C,  then 

Av  =  vQ  V  .  (17) 

C 

Since  E=hv,  this  can  be  considered  a  shift,  S,  in  photon  enemy, 
given  by 

S  =  v  E  (18) 

~c  y 

The  effect  of  this  energy  shift  on  the  Mossbauer  cross 
section  can  be  closely  approximated  in  the  case  R  «  2k  0^, 

S<kOd,  by 
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<r  =  1  <r0  f  f'  (  r2  )  (19) 

2  (  s2+  r2  ) 

(Ref  30:208).  By  examining  this  expression  it  can  be  seen  that  any 
energy  shift  which  is  of  the  order  of  a  linewidth  will  significantly 
affect  absorption.  Consider,  for  example,  the  relative  velocity  re¬ 
quired  to  shift  the  energy  by  the  value  r,  one  linewidth,  in  an 
experiment  involving  Fe^.  By  (18), 

v  =  r  C  =  C  =  10 cm/sec. 

Ev  Q 

67 

If  Zn  is  used,  the  corresponding  velocity  will  be  of  the  order  of 
-5 

10  cm/sec. 

Acceleration.  The  experiment  by  Pound  and  Rebka  (Ref  19) 
verified  a  simple  approach  to  the  effect  of  acceleration  (or  equally, 
gravity)  on  the  Mossbauer  Effect.  A  gamma  photon  can  be  consider¬ 
ed  to  have  an  apparent  mass,  given  by  Einstein's  relationship 

9 

E  =  MC“.  In  a  gravity  field  this  mass  is  acted  upon  by  the  force, 

weight;  and  movement  of  weight  through  the  gravity  field  causes  an 

energy  shift  in  the  photon. 

S  =  Mg  (h2  -  hj)  =  E  g  (h2  -  h^  (20) 

C2 

As  before,  S  is  the  energy  shift;  g  is  the  gravitational  acceleration; 
and  h2  and  hj  are  the  initial  and  final  height  of  the  photon  above  the 
ground.  For  a  nominal  height,  100  meters,  S  is  only  0.03r.  There¬ 
fore,  unless  large  distances  or  large  accelerations  are  involved,  the 
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effect  of  acceleration  is  small  compared  to  the  effect  of  relative 
velocity. 

Temperature.  In  addition  to  the  effect  of  the  ambient  tem¬ 
perature  on  the  fraction  of  zero  phonon  events  given  in  (7),  there 
is  a  second  effect.  This  effect  is  an  energy  shift  in  the  gamma 
photons  due  to  a  difference  in  temperature  between  the  source  and 
the  absorber.  The  energy  shift  arises  from  the  second-order 
Doppler  effect,  which  can  be  recognized  in  (16)  as  follows: 

A  V  =  v0  V  (C  -  v)_1  (16) 

Expanding  (16)  in  a  binomial  series, 

A  w  =  y  (  v  +  v2  — -  )  .  (21) 

(  C~  ) 

The  shift  in  energy  arises  from  the  second  term  in  the 
parentheses.  If  the  source  and  absorber  are  at  different  temperatures, 
then  the  squares  of  the  velocities  of  the  lattice  vibrations  will  be  dif¬ 
ferent.  The  effect  is  measured  by  the  expression 

S  -  C,  AT  (22) 

Ev  2 A  C2 

where  S  is  the  energy  shift,  Cj  is  the  specific  heat  of  the  lattice, 

AT  is  the  temperature  difference  between  source  and  absorber,  and 
A  is  the  gram  atomic  weight  of  the  isotope  (Ref  22:275).  For  Fe  , 

S  is  only  0.005  r  for  a  temperature  difference  of  1°K.  Therefore, 
normal  methods  for  keeping  the  source  and  absorber  at  the  same 
temperature  should  be  adequate  for  any  but  the  most  precise  cxperi- 
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merits  with  Fe^. 

The  question  arises,  why  do  not  lattice  vibrations  induce  a 
first  order  Doppler  shift?  These  vibrations  produce  maximum  velo¬ 
cities  of  the  order  10^  cm/sec  (Ref  15:18),  and,  as  previously 
discussed,  the  Mossbauer  Effect  is  sensitive  to  relative  velocities 
of  the  order  10  cm/sec.  The  answer  is  found  in  the  fact  that 
the  period  of  the  lattice  vibrations  is  of  the  order  10“13  sec  (Ref 
5:57).  The  effects  due  to  the  vector,  velocity,  tend  to  cancel  out 
over  the  "long"  lifetime  of  the  nuclear  excited  state  which  is  of  the 
order  10“^  seconds  in  Fe^?.  On  the  other  hand,  the  effects  due  to 
the  scalar  term,  v^,  are  additive. 

Chemical  Shifts.  If  the  source  and  absorber  differ  slightly 
in  chemical  or  isotopic  *  composition,  this  will  result  in  a  difference 
in  the  nuclear  transition  energy  involved  in  the  Mossbauer  Effect. 

As  a  consequence,  no  source -absorber  pair  will  exhibit  maximum 
Mossbauer  absorption  at  exactly  zero  relative  velocity,  zero  temp¬ 
erature  difference,  and  zero  acceleration.  In  addition,  no  two 
source -absorber  pairs  will  have  their  absorption  maxima  displaced  by 
exactly  the  same  amount  on  an  energy  scale.  However,  if  the  source 
and  absorber  are  carefully  prepared,  the  chemical  shift  can  be  mini¬ 
mized  to  within  a  small  fraction  (10-^)  of  r  (Ref  19:339). 

*  The  term  "chemical  shift"  is  used  in  the  literature  to  denote 
both  isotopic  and  chemical  effects  (Ref  8:53). 
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Magnetic  Effects.  Application  of  a  magnetic  field  to  either 
the  source  or  absorber  causes  splitting  of  the  energy  levels  in 
the  nuclear  excited  state.  The  net  result  is  that  there  will  be 
several  peaks  in  the  plot  of  the  Mossbauer  cross  section  versus 
energy,  instead  of  the  one  distinct  peak  shown  in  Figure  3.  For 
Fe^  the  magnitude  of  the  magnetic  field  necessary  to  produce 
noticeable  splitting  is  of  the  order  10^  gauss  (Ref  8:52). 
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III.  Limitations  on  the  Application  of  the  Mossbauer  Effect 

In  the  previous  Chapter  the  basic  theory  of  the  Mossbauer 
Effect  was  described,  and  relationships  showing  the  effects  of 
physical  phenomena  were  presented.  In  this  Chapter  an  analysis  will 
be  made  to  determine  the  order  of  accuracy  that  may  be  achieved 
in  measuring  physical  phenomena  by  use  of  the  Mossbauer  Effect. 

This  analysis  will  demonstrate  why  the  measurement  of  the  relative 
velocity  shift  is  the  most  promising  mechanism  for  practical  applica¬ 
tions,  and  why  the  choice  of  the  isotope  Fe^  is  most  clearly  indicated 
for  such  applications.  A  brief  discussion  of  other  important  limitations 
is  included. 

Counting 

Chapter  II  contained  an  explanation  of  how  a  change  in  the 
Mossbauer  cross  section  could  be  detected  by  counting  the  number 
of  photons  passed  through  an  absorber.  The  questions  that  must  be 
answered  are  how  accurately  can  the  change  in  cross  section  be 
measured,  and  how  accurately  can  the  change  in  cross  section  be 
related  to  the  physical  phenomenon  causing  the  change.  In  order  to 
answer  the  first  question,  it  is  necessary  to  describe  the  best  method 
of  counting  gamma  rays. 

Gamma  rays  are  most  efficiently  counted  by  a  scintillation 
detector,  shown  schematically  in  Figure  5.  The  photon  is  absorbed 
in  the  scintillator  by  Compton  scattering  or  the  photoelectric  process. 
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Schematic  Diagram  of  a  Scintillation  Detector  (Ref  23:163) 
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The  light  flashes  from  subsequent  atomic  transitions  are  converted 
to  electronic  pulses  by  the  photo- multiplier  tube,  which  also  amplifies 
these  pulses  many  times  (Ref  23:163).  The  pulses  can  then  be 
further  amplified,  measured,  and  counted  by  electronic  devices.  The 
use  of  a  discriminator  allows  elimination  of  all  pulses  except  those 
resulting  from  incident  radiation  in  a  particular  energy  range.  In 
addition  to  high  efficiency  (number  of  photons  counted  divided  by  the 
number  of  photons  incident  on  the  detector),  the  use  of  the  scintil¬ 
lation  detector  allows  the  highest  possible  counting  rate  (Fief  23:162). 

The  efficiency  of  the  detector  can  be  estimated  from  the 
number  of  photons  absorbed  by  the  scintillator.  A  sodium  iodide, 
thallium  activated  -  Nal(Tl)  -  crystal  less  than  0.1  cm  thick  will 
absorb  nearly  all  the  incident  photons  which  are  in  the  energy  range 
of  interest  (10  -  lOOkev)  for  Mossbauer  experiments  (Ref  23:190,  191). 
The  maximum  counting  rate  is  limite  i  by  the  resolving  time  of  the 
detector,  which  is  the  minimum  time  which  must  elapse  between 
successive  events  if  they  are  to  be  counted  as  two  events.  The 
resolving  time  is  determined  by  the  recovery  time  of  the  electronic 
scalar  and  by  the  decay  time  of  the  scintillator.  A  low  value  for 
resolving  time  is  10"^  seconds  (Ref  23:198).  This  would  allow  a 
counting  rate  of  10®  per  second  if  a  10%  correction  is  made  for 
counts  lost  during  the  dead  time  of  the  detector  (Ref  23:127). 

The  maximum  counting  rate  determines  a  basic  limitation  in 
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the  accuracy  of  a  Mossbauer  experiment.  This  results  from  the 
theory  of  statistics  of  detection  systems.  If  n  is  the  number  of 
counts,  then  the  standard  deviation,  an>  is  given  by 

crn  - 100  .  %  (23) 

(n)l/2 

(Ref  23:59).  If  the  standard  deviation  is  taken  as  an  approximate 
criterion  of  accuracy,  and  if  the  maximum  n  is  10°,  then  the  mini¬ 
mum  detectable  change  in  the  counting  rate  is  0.1%.  (For  the  purposes 
of  this  exercise,  it  is  assumed  that  it  is  desirable  to  obtain  a  data 
point  in  one  second;  therefore,  if  the  maximum  counting  rate  is  10®/ 
second,  n  =  10°.  ) 

Using  this  value,  it  is  possible  to  estimate  the  minimum  change 
in  the  Mossbauer  cross  section  that  can  be  detected.  A  parameter, 

F,  is  defined  as  the  ratio  of  the  flux,  I(x),  of  photons  that  reach 
the  detector  to  the  flux,  Io,  of  photons  incident  upon  the  absorber. 

F  =  I(x)  =  exp  (-  No  a  x)  (24) 

Io 

(Ref  13:344).  No  is  the  number  of  atoms/cmS  in  the  absorber  that 
are  the  proper  isotope  for  Mossbauer  absorption,  and  x  is  the  ab¬ 
sorber  thickness.  The  equation  assumes  that  the  only  absorption 
that  takes  place  is  due  to  the  Mossbauer  Effect,  that  the  flux  due 
to  re  radiation  is  negligible,  and  that  all  photons  which  pass  through 
the  absorber  reach  the  detector.  These  conditions  can  be  approximated 
by  a  suitable  geometry,  a  thin  absorber,  and  an  isotope  such  as  Fe^ 
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for  which  the  Mossbauer  cross  section  is  much  larger  than  the 
competing  cross  sections.  If  (24)  is  differentiated  with  respect  to 
a,  and  if  AF=dF,  then 

AF  =  -  No  x  Ag  exp  (  -No  ct  x)  (25) 

In  order  to  get  as  large  a  change  in  counting  rate  as  possible,  it 

is  desirable  to  maximize  A  F.  This  can  be  done  with  respect  to 

x,  the  absorber  thickness.  The  result  is  an  optimum  thickness 

x  (optimum)  -  1  (26) 

No  CT 

(An  optimum  iron  absorber  containing  a  natural  abundance  -  2%  -  of 
Fe^  would  be  about  0.05  cm  thick.)  If  (26)  is  substituted  in  (25), 

A  F  =  -  .368  Aa  (27) 

CT 

The  conclusion  can  now  be  realized  that  if,  by  previous  argument, 

the  maximum  change  in  counting  rate  that  can  be  detected  is  0. 1%, 

then  the  fractional  change  in  cr  that  can  be  detected  is  about  0.  3%. 

It  is  now  a  relatively  simple  matter  to  apply  these  results 

to  a  determination  of  the  minimum  detectable  change  in  the  energy 

of  a  physical  phenomenon.  The  necessary  equation  is  repeated: 

a-  1  CT0ff'  (_r2  )  (19) 

2  (  +  r2  ) 

It  is  convenient  to  let  S  =qr.  If  (19)  is  differentiated  with  respect 
to  q,  and  Act  =  dcr , 

A  CT=  1  CTo  f  f'  A  q  r  _ 2q  7  (28) 

2  ~  (1+  q2)2  ~ 
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By  substituting  (19)  and  (28)  in  (27), 

AF  =  Aq  /~  -  2q  ~]  (29) 

Tm?) 

(29)  can  be  maximized  with  respect  to  q,  giving  q  =  +  1.  Then, 
using  q  =  -1, 

AF  =  Aq  (30) 

This  means  that  the  minimum  detectable  change  in  the  energy  of  a 
physical  phenomenon,  when  expressed  in  units  of  r,  is  of  the  same 
order  as  the  minimum  detectable  change  in  the  counting  rate.  That 
is,  Aq  =  0.1%,  or  AS  =  10"^  r. 

Attention  should  be  called  to  the  fact  that  by  letting  q  =  -1, 
we  imply  that  the  emitted  photons  are  biased  in  energy  by  one  r 
away  from  the  central  maximum  of  the  absorption  spectrum,  i.e. 
on  that  part  of  the  o  vs  E  curve  where  the  slope  is  greatest.  An 
excellent  example  of  this  technique  was  provided  by  the  Pound  and 
Rebka  experiment  (Ref  19),  in  which  the  emitted  photons  were  shift¬ 
ed  in  energy  by  a  constant  velocity  drive  in  order  to  allow  the  small 
AS  due  to  the  gravity  effect  to  result  in  the  largest  possible  change 
in  counting  rate. 

One  further  consideration  is  whether  or  not  a  scattering  ex¬ 
periment  offers  an  advantage  in  the  matter  of  counting  limitations. 

In  an  ideal  scattering  experiment,  the  detectors  are  completely 
shielded  from  the  source,  but  detect  all  reradiation  from  the  absorber. 


26 


GA/Phys/62-7 

Let  the  re  radiated  flux  be  I',  where 

r  Io  -  I(x)  =  Io  r  1  -  exp  (-Noa  x)  7  (31) 

The  quantity  F’  is  defined, 

F'=  1=1-  exp  (-No a  x)  (24a) 

To" 

If  (24a)  is  differentiated  with  respect  to  a ,  an  equation  of  the  same 
form  as  (25)  is  obtained.  Therefore,  the  analysis  for  the  scatter¬ 
ing  experiment  is  the  same  as  for  the  absorption  experiment. 

However,  the  scattering  experiment  poses  certain  practical  difficulties 
since  reradiation  is  isotropic. 

Line  width 

In  the  section  just  presented,  it  was  demonstrated  that  the 
statistics  of  counting  limited  the  minimum  detectable  change  in 
energy,  AS,  of  a  physical  phenomenon  to  a  certain  fraction  of  the 
Mossbauer  linewidth,  r.  The  example  given  was  for  one  second  of 
counting  at  the  high  counting  rate  of  10®/  second,  in  which  case 
A S=  10'®  r.  The  next  questions  of  importance  are,  how  small  is 
r,  and  what  does  this  mean  in  terms  of  the  physical  phenomena 
pertinent  to  the  investigation. 

R7 

As  discussed  in  Chapter  II,  Zm  exhibits  the  smallest  line- 
width.  The  difficulty  in  using  Zn^'7  is  that  the  Mossbauer  cross  section 
is  so  small  as  to  prohibit  the  use  of  the  optimum  absorber  thickness 
given  by  (26).  If  one  uses  the  values  from  \ppendix  B,  calculation 
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will  show  that  even  for  an  absorber  enriched  to  100%  of  the  proper 
isotope,  x  (optimum)  for  Zn®7  is  about  13  cm.  In  order  to  pene¬ 
trate  an  absorber  this  thick  with  photons  of  the  100  kev  range,  and 
still  have  10®  photons/sec  for  counting  purposes,  will  require  an 
incident  flux  of  10^6  photons/second.  (The  calculation  is  made  using 
a  photoelectric  cross  section  only.)  To  provide  such  an  incident 
flux  would  require  a  source  strength  of  at  least  10^  curies,  which 
is  beyond  the  realm  of  reason  for  safety  considerations  alone.  An 
absorber  only  1  cm  thick  will  mean  in  terms  of  the  previous  analysis 

that  the  AS  which  may  be  detected  is  just  equal  to  r.  This  means 

67 

that  the  best  that  can  be  done  with  Zn  in  a  short  counting  time  is 
to  tell  whether  or  not  resonant  absorption  is  present.  As  a  con¬ 
sequence  of  this,  and  the  additional  difficulties  that  the  use  of  Zn67 
implies,  (very  low  temperatures  and  equipment  free  from  vibrations 
of  10 cm/sec)  it  is  not  deemed  advisable  tc  attempt  a  practical 
application  with  this  isotope. 

In  contrast  to  Zn67,  Fe57  has  a  large  enough  Mossbauer 
cross  section  to  permit  the  use  of  very  thin  absorbers.  Since  Fe^7 
also  exhibits  the  second  smallest  r  of  any  Mossbauer  isotope,  it  is 
appropriate  to  use  the  Fe^7  line  width  as  the  basis  for  calculations 
of  what  may  possibly  be  done  in  the  measurement  of  physical  phenom¬ 
ena. 
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It  was  calculated  in  Chapter  II  that  a  relative  velocity  of 
10‘2  cm/sec  will  shift  the  Fe^  absorption  spectrum  by  one  line- 

O 

width.  If  a  shift  of  10“°  r  is  measurable,  then  a  relative  velocity 
of  10 "5  cm/sec  may  be  detected. 

The  question  of  whether  or  not  the  Mossbauer  Effect  can  be 
used  directly  to  measure  acceleration,  A,  may  be  answered  by  re¬ 
placing  g  by  A  in  (20).  For  a  nominal  distance,  h2  -  hj_,  equal  to 
1  meter,  and  basing  calculation  on  Fe^?, 

S=3xlO'5rxA  (32) 

where  A  is  in  meters/sec^.  Again  using  the  minimum  detectable 

o 

AS  =  10  r,  the  measurable  change  in  acceleration  is  30  meters/ 
sec2,  or  about  3g!  (In  order  to  reduce  this  to  0.003g  would  require 
a  source -absorber  separation  of  1000  meters  or  a  counting  time  of 
10®  seconds.) 

Similar  calculations  can  be  carried  out  for  other  physical  phe¬ 
nomena  by  using  the  relationships  from  Chapter  II.  The  effects  of 
temperature  difference,  chemical  shifts,  and  magnetic  effects  are 
small  compared  to  the  effect  of  relative  velocity.  Therefore,  re¬ 
lative  velocity  is  chosen  as  the  mechanism  for  applying  the  Mossbauer 
Effect  to  accelerometers.  Techniques  to  achieve  this  end  are  de¬ 
veloped  in  Chapter  IV  and  Chapter  V. 

Other  Limitations 

Problem  areas  in  the  application  of  the  Mossbauer  Effect  can 
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be  recognized  from  the  effects  of  physical  phenomena  previously 
described.  Temperature  control  must  be  employed,  very  strong 
magnetic  fields  must  be  avoided,  and  care  must  be  used  to  prepare 
the  source  and  absorber  in  order  to  minimize  chemical  shifts.  There 
is  no  serious  shielding  problem  in  the  use  of  Fe^.  A  few  milli¬ 
meters  of  material  will  reduce  the  gamma  flux  to  an  acceptable 
level  unless  an  extremely  strong  source  is  employed.  One  of  the 
most  serious  problems  is  that  of  vibration.  Small  vibrations  result¬ 
ing  in  source -absorber  relative  velocities  must  be  avoided.  This 
could  be  a  near  impossible  obstacle  for  a  vehicular  operation. 

The  best  statement  that  can  be  made  with  regard  to  most  of 
these  limitations  is  that  they  have  been  overcome  in  the  precise 
scientific  experiments  already  mentioned.  One  would  expect  that  as 
these  experiments  continue,  techniques  will  be  perfected  and  knowledge 
of  the  basic  physics  will  increase.  With  regard  to  vehicular  vibrations, 
rigid  source -absorber  coupling  may  minimize  this  difficulty.  In 
space  flights  using  low  thrust  propulsion  devices  the  problem  should 
be  much  less  severe  than  for  flights  in  the  atmosphere.  In  the  re¬ 
mainder  of  this  study,  the  assumption  will  be  made  that  these  problems 
can  be  overcome  by  diligent  research  and  engineering.  The  basic 
limitations  of  counting  and  the  line  width  of  Fe^  will  be  used  to  define 
the  possibilities  of  the  Mossbauer  Effect.  If  future  discoveries  provide  a 
better  isotope,  the  results  of  the  study  may  be  easily  extrapolated. 
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IV.  Techniques  for  Measurement  o£  Velocity  and  Position 

As  discussed  in  Chapter  III,  it  is  impractical  to  measure 
acceleration  by  the  direct  use  of  the  Mossbauer  Effect.  On  the 
other  hand,  the  Mossbauer  Effect  is  very  sensitive  to  the  first- 
order  Doppler  shift  produced  by  a  small  relative  velocity  between 
the  source  and  the  absorber.  In  order  to  apply  this  phenomenon 
to  accelerometers,  it  is  necessary  to  find  techniques  by  which  the 
extreme  sensitivity  of  the  Mossbauer  Effect  to  small  velocities  can 
be  used  as  an  instrument  sensor  or  pickoff  device.  Two  such 
techniques  are  presented.  One  technique  allows  measurement  of 
velocity.  Another  provides  for  measurement  of  angular  position. 
Chapter  V  contains  an  explanation  of  how  these  techniques  can  be 
applied  to  accelerometers. 

Measurement  of  Velocity 

Direct  measurement  of  velocity  by  the  use  of  the  Mossbauer 
Effect  is  impractical  in  many  instances  because  only  small  velocities 
can  be  measured.  Two  systems  are  presented  which  allow  measure¬ 
ment  of  velocities  of  any  magnitude.  The  first  of  these  systems 
will  be  referred  to  as  the  phase  difference  method.  The  phase  dif¬ 
ference  method  is  applicable  to  measurement  of  angular  velocity. 

Phase  Difference  Method.  As  shown  in  Figure  6a,  emitter 
atoms  are  placed  at  discrete  points  1  and  2  on  the  rim  of  a  rotating 
wheel.  When  the  wheel  is  rotating  through  the  interval  -  0Q  to  +  0o, 
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Phase  Difference  Method 
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the  emitted  photons  passing  in  a  direct  line  to  the  absorber  have 
a  velocity  component,  ve,  along  this  line  equal  to  Ve  sin  0=Ve  sin 
cut.  Ve  is  the  tangential  velocity  of  emitter  atoms,  cu  is  the  angul¬ 
ar  velocity  of  the  wheel.  ve  is  plotted  for  both  emitter  1  and 
emitter  2  in  Figure  6b.  Motion  is  induced  in  the  absorber  through 
the  vibration  of  an  electromagnetic  core.  In  practice,  the  absorber 
foil  may  be  mounted  in  a  very  narrow  slit  in  this  core.  The  width 
of  the  slit  will  subtend  and  define  the  arc  2  0o.  Only  gamma  rays 
following  the  direct  line  normal  to  the  absorber  will  be  significant 
in  the  counting,  since  the  photons  striking  the  core  or  surrounding 
apparatus  will  be  absorbed.  The  relative  velocity  of  the  absorber 
is  assumed  to  be  sinusoidal. 

Let  the  absorber  velocity,  va,  be  Va  sin  wt  where  Va  is 
the  amplitude  of  the  absorber  velocity  and  w  is  the  angular  fre¬ 
quency.  We  shall  require  initially  that  relative  velocity,  vr, 
between  the  source  and  absorber  be  as  small  as  possible  in  the 
interval  -  0O  <  at  <+0o.  Then,  for  this  interval,  let 

vr=  0  =  Ve  sin  cut  -  Va  sin  wt  (33) 

or,  if  0o  is  small,  Ve  c u  =  Va  w. 

One  possible  solution  is:  Ve  =  Va,  <u  =  w. 

However,  since  0o  is  small,  this  will  provide  a  minimum  of  counting 

time  during  one  complete  revolution  of  the  wheel.  When  a  second 

emitter  is  used,  it  is  necessarv  to  let  Ve  =  2  Va,  tu=  w  .  Thus,  the 

T 
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counting  time  available  is  doubled.  As  illustrated  in  Figure  6b, 

the  absorber  velocity  is  nearly  coincident  with  the  emitter  velocity 

every  half-cycle  of  revolution.  Subject  to  limitations  which  will 

be  presently  examined,  we  may  further  increase  counting  time  by 

placing  more  emitters  on  the  wheel  rim  and  by  requiring,  in  general 

Ve  =  n  Va,  o>=w  ,  where  n  is  the  number  of  emitters.  If  Fe^  is 
n 

used  as  the  emitter  isotope,  an  appreciable  Mossbauer  absorption 
will  take  place  as  long  as  vr  is  less  than  10"2  cm/sec.  The  out¬ 
put  of  the  instrument  can  be  the  electronic  frequency,  w,  of  the 
excitation  voltage  used  to  vibrate  the  electromagnetic  core.  The 
fact  that  w  is  larger  than  u>  presents  an  advantage  since  in  measur¬ 
ing  cyclic  events  we  may  expect  an  accuracy  of  at  least  +  one  cycle. 
More  generally,  if  we  can  measure  a  given  part  of  one  cycle,  then 
the  more  cycles  available,  the  more  accurate  will  be  our  measurement. 

In  order  to  establish  the  sensitivity  of  this  technique  it  is 
assumed  that  at  some  time,  -tQ,  the  absorber  velocity  is  in  phase 
with  the  emitter  velocity.  Subsequently  a  small  error  in  frequency 
of  the  absorber,  aw,  results  in  a  phase  difference  at  time  after  -tQ. 
This  is  illustrated  in  Figure  6c.  The  resultant  relative  velocity  may 
be  expressed  as  follows  for  the  interval  -o0  <0  <+0o. 

vr  =  ^e  sin  “  ^a  s*n  (n  +  ^  (34) 

where  A  is  the  phase  difference. 

6  =  n  A  co  (t  +  t0)  (35) 
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9 

If  e0  is  very  small,  say  10  radians,  and  to  is  some  nomin¬ 
al  value,  say  10  radians  per  second,  then  t  is  only  10"^  seconds. 

On  the  other  hand,  t  should  be  on  the  order  of  one  second  if 
enough  counts  are  to  be  obtained  to  establish  whether  there  is  a 
relative  velocity  between  emitter  and  absorber.  Therefore,  the 
approximation  is  made,  for  t  «t0, 

<t>  =n  A  toto  (36) 

Substituting  (36)  in  (34), 

vr=Ve  sin  tot  -  Va  sin  n  (tot+Atot0)  (37) 

or, 

vr=  Ve  sin  tot  -  Va  sin  n  tot  cos  n  Awt0 

-  Va  sin  n  AtotQ  cos  n  tot  (38) 

If  Va  =  ,  and  if  all  the  arguments  are  small  angles, 

n 

vr=-  Ve  Ato  tQ  (39) 

The  smallest  increment  of  angular  velocity  that  can  be  measured  is 
therefore, 

A  to  =  vf  (40) 

Veto 

where  vr  is  the  relative  velocity  that  can  be  detected  by  the  Moss- 
bauer  Effect.  If  vr  =  10'^  cm  /sec,  t0  =  1  sec,  and  Ve  =  Rto  = 

n 

100  cm /sec,  then  an  increment  of  10  radians/sec  will  be  measur¬ 
able.  An  interesting  feature  of  this  technique  is  that  the  accuracy 
increases  for  large  values  of  to. 
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The  above  description  was  based  on  ideal  conditions  of  opera¬ 
tion.  It  is  pertinent  to  examine  some  of  these  conditions  in  more 
detail  in  order  to  establish  the  limitations  of  the  technique. 

1.  It  was  assumed  that  the  amplitude  of  the  absorber 
velocity  is  set  at  precisely  the  right  value.  The  relative  velocity 
due  to  amplitude  deviation,  A  Ve,  is 

vr  =Ve  sin  cut  -  (  Ve  +  A  Ve  )  sin  n  cut  (41) 

(  n  ) 

This  is  at  a  maximum  during  the  interval  of  counting  when  cut  =  G0. 


(max)  vr  =  -  A  Ve  0o 

The  relative  velocity  due  to  amplitude  error  should  be  less  than 

the  relative  velocity  due  to  phase  difference  in  order  for  this  technique 

to  be  effective  .  Therefore,  it  is  required  that 

A  Ve  0o  <  Ve  A  cutG  (42) 


If  t0  =  1  sec 


A  Ve  <,  A  cu  (43) 

^e  °o 

If  A  cu  is  to  be  measured  to  an  accuracy  of  10-3  radians  /sec  and 
Qq  =  10-3  radians,  then  the  amplitude  of  the  absorber  velocity  must 
be  held  accurately  to  less  than  one  part  in  100.  The  relationship 
points  out  the  advantage  of  making  0O  as  small  as  possible. 

2.  Another  reason  for  making  e0  small  is  to  allow  the  use 
of  as  many  emitters  as  possible,  thus  making  the  frequency  multi¬ 
plication  factor,  n,  as  large  as  possible.  The  question  raised  is,  if 
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or 


vr  =  Ve  sin  ut  -  Ve  sin  nut 

*17 


(44) 


vr  (max)  =  Ve  sin  90  -  sin  n  G0  (45) 

n 

how  small  must  0O  be  to  allow  a  given  number  of  emitters,  n,  with¬ 
out  causing  the  relative  velocity  from  (45)  to  be  larger  than  the  re¬ 
lative  velocity  induced  by  the  phase  error?  This  problem  may  be 
analyzed  by  expanding  both  terms  on  the  right  hand  side  of  (45)  in 
a  series,  and  combining  the  series  term  by  term,  yielding: 

v r  (max)  =  Ve  /"  (n2  -  1)  0O3  -  (n4  -  1)  o05+ (n6  -  l)  0o7---7(46) 

^ - ^ -  - 


It  is  required  that  the  relative  velocity  due  to  phase  error  (39)  be 
larger  than  the  vr  (max)  given  by  (46),  or 

Ao>t0  >  /~  (n2  -  l)0o3~— 7  (47) 

- 3TT^ 


Again  it  is  assumed  tQ  =  1  second,  allowing  computation  of  approxi¬ 
mate  values.  A  cu  must  be  numerically  greater  than  the  figures 
represented  in  the  elements  of  the  Table  I.  A  typical  result  of  the 
computation  is  that  if  a  large  number  of  emitters  are  to  be  used, 
then  oQ  must  be  on  the  order  of  10"^  radians  in  order  that  a  fre¬ 
quency  error  on  the  order  of  10-^  radians  /  second  can  be  measured. 

3.  There  is  a  minimum  below  which  0O  cannot  be  reduced. 
This  limit  is  established  by  the  physical  dimensions  of  the  wheel.  The 
arc  subtended  on  the  wheel  by  the  emitter  must  be  less  than  the  arc 
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Table  I 

Limiting  Increments  of  Angular  Velocity 

0O  n  =  2  n  =  8  n  =  16  n  =  100 

.  1  rd  .5x10"^  1.0  x  10-2  -  - 

.01  rd  .  5  x  10-6  1.05  xlO”5  4.25  x  10-5  1.6  x  10-3 

.  001  rd  .  5  x  lO'9  1.  05  x  10“8  4.  25  x  10“8  1. 67  x  10-6 

subtended  by  0o  in  order  for  the  foregoing  calculations  to  have  mean¬ 
ing  (a  point  source  was  assumed).  In  the  examples  previously 
discussed,  R=10cm.  For  9o=10-8  radians,  the  arc  subtended  is 
10  cm.  This  would  indicate  that  the  arc  subtended  by  the  emitter 
should  be  10“^  or  10'5  cm,  or  on  the  order  of  the  wavelength  of 
visible  light.  If  a  number  of  emitters  are  used,  then  the  spacing 
between  emitters  should  be  accurate  to  10“^  or  10“®  cms.  This 
will  require  precise  and  expensive  manufacturing  techniques.  How¬ 
ever,  the  requirements  should  not  be  declared  impossible  in  view 
of  progress  in  solid  state  physics.  For  instance,  superconductivity 
experiments  have  been  carried  out  with  films  10"  ^  cm  in  thickness 
(Ref  10:459). 

4.  It  is  evident  from  previous  discussion  that  the  length 
of  time  available  for  measurement  is  a  critical  limitation  of  this 
technique.  Time  is  necessary  to  obtain  statistical  counting.  Time 
(tQ)  is  necessary  to  allow  a  small  frequency  error  to  result  in  a 
phase  difference  between  absorber  and  emitter.  Time  is  necessary 
to  count  a  sufficient  number  of  cycles  of  the  coil  excitation  voltage, 
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in  order  to  read  the  angular  velocity  to  the  desired  accuracy. 

Finally,  since  the  phase,  frequency,  and  amplitude  of  the  absorber 
must  be  accurately  adjusted,  there  may  be  a  significant  time  neces¬ 
sary  for  even  an  electronic  feedback  system  to  establish  "zero" 
relative  velocity  between  source  and  absorber. 

The  four  limitations  just  described  establish  guidelines  for 
the  use  of  the  phase  difference  method  for  applying  the  Mossbauer 
Effect  to  measurement  of  angular  velocity.  A  second  velocity 
measurement  technique  will  now  be  examined.  This  technique  is 
referred  to  as  the  amplitude  difference  method.  There  are  two 
variations  of  the  system.  The  first  variation  employs  discrete 
emitters  and  is  applicable  to  measurement  of  angular  velocity.  The 
second  variation  employs  continuous  emission  and  is  applicable  to 
measurement  of  both  angular  velocity  and  linear  velocity. 

Amplitude  Difference  Methods.  Figure  7a  is  a  schematic 
representation  of  the  amplitude  difference  method  using  discrete 
emitters.  Point  source  emitters  are  placed  at  positions  one  and 
two  on  the  rim  of  the  wheel  which  is  rotating  through  an  angle  e, 
measured  from  a  line  parallel  to  the  plane  of  the  absorber.  The 
absorber  is  mounted  in  a  slit  in  an  electromagnetic  core  similar  to 
the  device  used  in  the  phase  difference  method.  (In  practice,  relative 
velocity  of  the  absorber  may  be  provided  by  a  number  of  techniques; 
i.e.,  piezoelectric  effect,  small  wheel  driven  by  a  synchronous 
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motor,  tuning  fork,  etc.)  The  edge  of  the  shielding  below  the  ab¬ 
sorber  (bottom  edge  of  the  slit)  defines  one  point  on  the  arc  0O. 

The  other  point  is  defined  by  the  tangent  to  the  top  of  the  wheel. 
fl0  is  the  interval  during  which  photons  from  the  emitter  are  counted 
by  the  detector.  If  0o  is  very  small,  only  photons  that  pass  on  the 
line  perpendicular  to  the  absorber  will  be  counted.  In  the  interval 
of  emission,  0<0<9o,  the  relative  velocity  of  the  emitter  is 

ve  =ve  cos  (48) 

where  Ve  =  Rco,  the  tangential  velocity  of  the  rim  of  the  wheel. 

The  velocity  of  the  absorber  is  sinusoidal,  given  by 

va  =  Va  cos  2wt  (49) 

for  the  two  emitters  shown.  In  general,  for  n  emitters,  va  = 

Va  cos  nwt.  For  case  I  plotted  in  Figure  7b,  there  is  a  relative 
velocity  between  emitter  and  absorber: 

vr  =Ve  cos  cut  -  Va  cos  2  cut  (50) 

The  minimum  vr  occurs  if  Ve  =  Va.  If  this  condition  is  fulfilled, 
the  minimum  relative  velocity  results  in  a  maximum  of  Mossbauer 
absorption,  and  the  angular  velocity  of  the  wheel  may  be  measured 
by  an  electronic  reading  of  the  frequency  or  amplitude  of  the  excita¬ 
tion  current  of  the  electromagnetic  core.  If  this  condition  is  not 
fulfilled,  there  will  be  a  larger  relative  velocity  difference  as  shown 
in  Figure  7b,  Case  II.  If  the  angles  are  small  in  (50), 

vr  ~  ve  “  va  =  (51) 
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Since  A  V  =  R  Aco 

Au)  =  vr  (52) 

TT 

For  vr  =  10"^  cm /sec,  R  =  10  cm,  an  angular  velocity  may  be 
measured  to  10"®  radians/second. 

Much  of  the  discussion  regarding  ihc  limitations  of  the  phase 
difference  method  is  also  applicable  to  the  amplitude  difference 
method.  The  major  difference  is  that  in  seeking  to  measure  the 
velocity  resulting  from  a  difference  in  amplitude  of  the  velocity 
cycles,  we  incur  an  obligation  to  accurately  control  the  amplitude 
of  the  absorber  velocity.  Further,  if  a  frequency  output  is  assum¬ 
ed  to  be  desirable,  it  is  necessary  to  accurately  relate  the  frequency 
of  the  absorber  vibration  to  the  maximum  absorber  velocity.  For 
a  vibrating  electromagnetic  core  this  relationship  is  established 
through  the  length,  1Q,  through  which  the  core  travels  in  one -half 
cycle;  i.e.,  if 


x  =  lo  sin  cot 

(53) 

dx  =  ln  oi  cos  co  t 
cTF 

(54) 

For  an  electromagnetic  core,  1Q,  may  be  a  variable  due  to  effects 
of  friction,  hysterisis,  electronic  noise,  and  other  causes.  On  the 
other  hand,  for  a  wheel  (  10  =  radius),  1Q  is  a  fixed  quantity.  There¬ 
fore,  for  the  amplitude  difference  method  it  may  be  desirable  to 
use  small  wheels  driven  by  a  synchronous  motor,  with  absorbers 
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mounted  as  shown  in  Figure  8.  The  output  quantity  now  becomes 
the  frequency  of  the  synchronous  motor.  If  n  emitters  are  used, 
then  there  is  a  frequency  multiplication  factor,  n,  between  the 
angular  frequency  of  the  emitter  wheel  and  the  frequency  of  the 
smaller  absorber  wheels.  The  question  may  be  raised,  why  not 
measure  the  velocity  of  the  emitter  wheel  directly  by  use  of  a  syn¬ 
chronous  motor  drive?  Other  than  the  frequency  multiplication 
factor  just  mentioned,  there  is  no  reason  why  this  cannot  be  done, 
if,  in  fact,  it  is  permitted  by  the  nature  of  the  physical  system 
involved.  There  are  certain  applications  wherein  it  is  necessary 
to  allow  a  wheel  to  rotate  as  freely  as  possible.  In  such  instances 
it  is  not  permissable  to  drive  the  emitter  wheel. 

The  number  of  emitters  that  may  be  employed  in  order  to 
obtain  a  frequency  multiplication  factor  is  limited  by  the  size  of 
the  angle  0o  and  the  tangential  velocity  of  the  emitter,  Ve.  If 


(max)  vr  =  Ve  cos  0o  -  Va  cos  n  0o  (55) 

and  if,  as  required,  Va  =  Ve,  then  by  expanding  the  cosine  terms 

in  series  and  subtracting  term  by  term 

(max)  vr  =  Ve  r  (n2  -  1)  0o2  -  (n4  -  1)  o04  +  (n6  -  1)  0o6  (56) 

21  31  B! 

It  is  imperative  that  (max)  vr  given  by  (56)  be  much  less  than  AV 

given  by  (51).  It  is  required  that 

A  Ve  >  T  (n2  -  1)  0o2  -  (n4  -  1)  0o4  +  (n^  -  1)  oQ6  - ~/  (57) 

Ve  21  31  K] 
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Figure  8 

Absorbers  Mounted  in  Rotating  Wheels 
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This  places  a  limit  on  the  accuracy  of  the  amplitude  difference 
method,  and  in  addition,  places  a  limit  on  the  magnitude  of  the 
velocity  which  may  be  measured.  There  was  no  such  restriction 
involved  in  the  phase  difference  method.  On  the  other  hand,  it  is 
possible  by  the  amplitude  difference  method  to  place  the  emitter 
atoms  along  the  entire  arc  0O.  Another  advantage  of  this  technique 
is  that  there  is  no  requirement  for  a  time  interval  in  which  a 
frequency  error  builds  into  a  phase  error. 

In  order  to  avoid  the  limitations  inherent  in  a  sinusoidally 
moving  absorber,  it  is  natural  to  question  whether  the  amplitude 
difference  method  can  be  employed  with  an  absorber  moving  with  a 
constant,  linear  velocity.  This  question  leads  to  a  consideration 
of  the  second  variation  of  the  amplitude  difference  method.  This 
variation  employs  continuous  emission  and  is  applicable  to  measure¬ 
ment  of  either  angular  or  linear  velocity.  As  applied  to  angular 
velocity,  the  method  is  illustrated  in  Figures  9a  and  9b.  Source 
atoms  are  placed  around  the  entire  rim  of  the  emitter  wheel.  Only 
those  photons  which  are  emitted  from  the  arc  of  the  rim  subtended 
by  the  angle  Oo,  and  of  these  only  those  photons  emitted  in  the 
direction  parallel  to  the  shielding  material  shown,  are  allowed  to  pass 
through  the  absorber  and  to  the  Detector.  These  photons  have  a 
relative  velocity  along  their  direction  of  propogation  given  by 

Vc  =  R  (jl>  COS  O  =  R  Cl)  (58) 
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Figure  9 


Amplitude  Difference  Method-Continuous  Emission 
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Relative  velocity  of  the  absorber  may  be  provided  by  either  of  the 
methods  shown  in  Figure  9.  In  Figure  9a  this  velocity  is  provided 
by  a  belt  drawn  by  small  wheels  of  radius,  r.  In  Figure  9b  the 
absorber  velocity  is  provided  by  the  rotation  of  a  wheel  tilted  at 
an  angle  p  to  the  direction  of  photon  propogation  .  Only  a  small 
section  of  the  wheel  is  used  as  an  absorber.  This  section  is  indicat¬ 
ed  by  the  darkened  portion  in  the  face  of  the  wheel  shown  in  Figure 
9b.  In  both  instances,  the  absorber  wheel,  or  wheels,  are  rotating 
at  angular  frequency  w.  This  frequency  may  be  controlled  by  a 
synchronous  motor.  The  output  of  the  system  may  be  the  frequency 
of  this  synchronous  motor.  In  both  cases,  the  velocity  of  the 
absorber  is  given  by 

va  -  r  w  sin  p  (59) 

The  relative  velocity  between  source  and  absorber  is  approximately 

vr-  R  u>-  r  w  sin  p  (60) 

By  this  technique  a  large  frequency  multiplication  factor  may  be 
obtained,  not  only  by  making  r  smaller  than  R,  but  by  making  p  a 
relatively  small  angle.  This  frequency  multiplication  factor,  n,  is 
given  by  the  condition  making  vr  zero.  Let  w  =  n  to,  then  by  (60) 

n  =  R  (61) 

r  sin  p 

The  accuracy  of  the  system  is  primarily  determined  by  the  size 
of  the  emitter  wheel.  A  deviation,  A  oj,  in  the  angular  frequency  of 
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the  emitter  wheel  will  result  in  a  relative  velocity, 

vr  =RAW  or  vr  =AVe  (62) 

A  limitation  on  the  system  accuracy  is  determined  by  the 

deviation  of  ve  due  to  the  size  of  0Q.  By  (58),  v0  will  vary  from 

R  u>  to  Rw  cos  0O.  It  should  be  noted  that  this  limitation  also 

applies  to  the  amplitude  difference  method  using  discrete  emitters, 

in  the  event  the  emitter  atoms  are  placed  along  the  entire  arc  0 

The  advantages  of  the  amplitude  difference  method  using  continuous 

emission  are  the  absence  of  any  inherent  error  due  to  frequency 

multiplication  and  the  absence  of  the  requirement  to  adjust  the  phase 

of  the  absorber  velocity.  These  advantages  should  make  it  a  simpler 

and  less  expensive  technique. 

In  summary,  for  the  phase  difference  method; 

A  u>=  vr  (40) 

Rajt0 

and  for  the  amplitude  difference  method; 

Av  =  Vj.  (52) 

R 

These  relationships  are  valid  only  within  the  limitations  discussed. 

A  plot  of  (40)  and  (52)  is  shown  in  Figure  10.  It  is  assumed  in 
this  plot  that  tQ  =  1  second,  R  =  10  cm. 

Measurement  of  Position 

A  technique  is  possible  for  using  the  Mossbauer  Effect  as  an 
accurate  null  sensor  for  angular  position.  As  shown  in  Figure  11 
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Figure  10 

Accuracy  of  Angular  Velocity  Measurement 


emitter  wheel 
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Measurement  of  Angular  Position  by  the  Mossbauer  Effect 
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emitter  atoms  are  placed  in  a  concentric  ring  of  radius,  r,  on 
the  face  of  the  emitter  wheel.  Shielding  is  employed  such  that 
only  the  atoms  in  the  ring  passing  an  aperture  in  the  shielding 
produce  photons  that  pass  through  the  absorber  and  to  the  detector. 
This  is  indicated  by  the  darkened  portion  on  the  face  of  the  wheel. 
The  photons  so  emitted  have  a  velocity  with  respect  to  a  stationary 
absorber,  given  by 

ve  =  rcu  sin  p  (63) 

If,  for  a  particular  application,  it  is  desired  to  maintain  the  wheel 
at  a  null  position  parallel  to  the  plane  of  the  absorber,  then  for  this 
application  we  may  take  p  as  a  small  angle.  The  accuracy  that  may 
be  obtained  is  given  by 

P  =  v§_  =  Yr_  (64) 

red  rid 

A  limitation  in  this  technique  is  the  variation  in  the  direction 
of  photon  propagation  due  to  a  need  for  a  finite  aperture  in  the 
shielding.  An  aperture  only  10"4  cm  wide  will  subtend  an  angle  of 
2"  of  arc  at  a  distance  of  10  cm.  Therefore,  the  smallest  angle 
that  could  be  detected  under  these  conditions  would  be  greater  than 
2"  of  arc.  An  aperture  of  10-^  cm  is  only  10,  000  Angstrom  Units 
and  is  equivalent  to  the  spacing  in  a  fine  diffraction  grating  for 
light  waves  (Ref  10:86). 
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V.  Application  to  Accelerometers 

In  Chapter  IV  techniques  were  presented  by  which  the  Moss- 
bauer  Effect  could  be  used  as  a  sensor  to  measure  the  velocity  or 
position  of  an  instrument  component.  In  this  Chapter  these  techni¬ 
ques  are  applied  to  accelerometers.  The  purpose  of  applying  the 
techniques  is  to  increase  the  resolution  or  to  decrease  the  thresh¬ 
old  sensitivity  of  accelerometers. 

Resolution  is  defined  as  the  reciprocal  of  the  least  fraction 
of  acceleration  which  can  be  measured.  That  is,  R=A/AA,  where 
R  is  resolution,  A  is  the  acceleration,  and  A  A  is  the  least  incre¬ 
ment  of  acceleration  which  can  be  measured.  Threshold  sensitivity 
is  the  smallest  detectable  unit  of  acceleration.  Threshold  sensitivi¬ 
ty  is  distinguished  from  the  term,  sensitivity,  which  is  the  ratio 
of  an  output  quantity  to  an  input  quantity  for  a  given  system  or 
component. 

In  addition  to  resolution  and  threshold  sensitivity,  time  re¬ 
sponse  or  frequency  response  is  an  important  performance  criterion. 
Time  response  indicates  the  highest  rate  of  change  of  acceleration 
that  can  be  measured.  Frequency  response  is  a  more  restricted 
statement  of  the  same  property.  Frequency  response  gives  the 
highest  cyclical  rate  of  acceleration  that  an  accelerometer  will 
read.  Instrumentation  using  the  Mossbauer  Effect  can  generally 
be  expected  to  consume  one  or  more  seconds  to  produce  a  single 
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data  point  because  of  the  need  for  statistical  counting.  Therefore, 
there  is  little  basis  for  improving  the  time  response  of  conven¬ 
tional  instruments,  which  may  measure  acceleration  of  40  or 
more  cycles/second  (Ref  27:87-88).  This  property  must  be  con¬ 
sidered,  however,  in  designing  any  accelerometer  and  the  time 
response  limitations  of  Mossbauer  techniques  should  be  appreciated. 

In  addition  to  resolution,  threshold  sensitivity,  and  time 
response,  many  other  terms  which  denote  the  performance  of  accel¬ 
erometers  are  found  in  the  literature.  These  terms  are  usually 
measures  of  limitations  rather  than  positive  indications  of  performance. 
Examples  are  non  linearity,  hysterisis,  friction,  mass  unbalance, 
and  temperature  effects. 

In  order  to  orient  the  study  that  follows,  it  is  worthwhile  to 
classify  accelerometers  in  several  major  categories.  Accelero¬ 
meters  from  each  category  will  then  be  analyzed  from  the  standpoint 
of  applications  using  the  Mossbauer  Effect. 

Classification  of  Accele rometers 

All  accelerometers  used  for  inertial  navigation  employ  a  mass, 
or  inertial  element,  which  is  sensitive  to  acceleration  through  the 
principle  of  Newtonian  mechanics,  force  equals  mass  times  accelera¬ 
tion.  The  response  of  the  inertial  element  may  be  measured  and 
quantitatively  related  to  acceleration  by  three  basic  methods.  These 
methods  serve  to  conveniently  classify  accelerometers: 
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1.  Simply  restrained  devices:  These  accelerometers 
use  the  elastic  property  of  a  material  or  other  medium  to  con¬ 
strain  the  motion  of  the  inertial  element.  Familiar  examples  are 
the  strain  guage  and  market  scale.  The  restraining  force  is 
proportional  to  an  elastic  constant,  K,  and  the  displacement,  x, 
of  the  inertial  element. 

2.  Force  balance  systems:  Constraint  is  generated 
through  a  servomechanism  which  attempts  to  maintain  the  inertial 
element  at  a  null  position.  The  constraining  device  is  usually  an 
electric  motor  or  other  form  of  electromagnetic  forcer.  The 
magnitude  of  the  constraint  is  proportional  to  thu  voltage  and/or 
current  supplying  the  forcer. 

3.  Unrestrained  devices:  The  inertial  element  is  allowed 
to  move  freely  under  the  influence  of  acceleration.  Output  is  a 
measurement  of  the  displacement  of  the  inertial  element  combined 
with  a  timing  device.  An  example  is  the  pendulum. 

Other  classifications  often  found  in  the  literature  include  the 
following: 

a.  Linear  or  pendulous  -  A  linear  accelerometer 
constrains  the  inertial  element  to  move  along  a  line.  A  pendulous 
accelerometer  constrains  the  sensitive  mass  to  describe  an  arc  of 
fixed  radius. 

b.  Singly  integrating  or  doubly  integrating  -  indicates 
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an  output  of  velocity  or  position,  respectively. 

c.  Solid,  liquid,  or  gas  -  may  indicate  the  phase  of 
the  inertial  element  and/or  the  constraining  device. 

d.  Other  -  such  adjectives  as  electromagnetic,  radio¬ 
active,  gyroscopic,  temperature-controlled,  aniso-elastic  .compensated, 
etc.  are  often  used  by  manufacturers  to  classify  accelerometers. 

These  adjectives  are  indicative  of  details  in  the  construction  of  ac¬ 
celerometers. 

Consideration  will  be  given  to  an  application  of  techniques  using 
the  Mossbauer  Effect  to  several  accelerometers  from  the  three  major 
classifications  described. 

Application  to  Simply  Restrained  Accelerometers 

Simply  restrained  devices  may  be  represented  schematically 
as  shown  in  Figure  12.  The  differential  equations  of  motion  are 
written  as  follows: 

For  Figure  12a,  a  linear  device, 

MA  =  Mx  +  Cx  +  Kx  (65) 

For  Figure  12b,  a  pendulous  device, 

MrA  =  IX  +  CX  +  KX  (66) 

The  steady-state  solution  of  (65)  is 

x  =  M  A  (67) 

IT 

The  resolution  of  a  simply  restrained  device  may  be  limited 
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a.  Linear 


I - 1 


A  -  input  acceleration,  cms/sec^" 

K  -  spring  constant,  dynes/cm 
C  -  viscous  friction  force,  dyne-sec/cm 
M  -  mass,  gms 
x  -  displacement  of  M,  cms. 


b.  Pendulous 


A  -  input  acceleration,  cms/sec^ 

K  -  elastic  constant,  dynes/radian 
C  -  viscous  friction  force,  dyne -sec/radian 
M  -  mass,  gms 

X  -  angular  displacement  of  M,  radians 
r  -  distance  from  pivot  to  center  of  M,  cms 
I  -  (not  shown)  -  moment  of  inertia  of  simple  pendulous 
element  (Mr^) 


Figure  12 

Schematic  Diagrams  of  Simply  Restrained  Accelerometers 


49a 


GA/Phys/62-7 

by  friction  or  by  uncertainty  in  the  value  of  K  caused  by  non 
linearity  and  hysterises  (Ref  32:Ch.  Ill,  pp.  72-81) .  An  example  is 
the  Whittaker  Z-125  Linear  Accelerometer  which  is  limited  to  end 
range  resolution  of  50  by  non  linearity  (Ref  27:84-85).  Another 
example  is  the  Modified  Seismic  System  by  Research,  Inc.,  which 
has  a  resolution  of  10  at  lg  input  (Ref  27:62-63).  In  addition  to 
being  limited  in  resolution,  simply  restrained  devices  are  limited 
in  threshold  sensitivity  by  static  friction  and  hysterisis.  It  will 
be  noted  in  (65)  that  it  is  necessary  to  have  a  significant  value  of 
C  in  order  to  damp  out  the  transient  solution.  It  is  difficult  to 
conceive  of  a  means  of  obtaining  high  C  without  accompanying 
static  friction. 

A  further  difficulty  arises  in  obtaining  an  output  from  a 
simply  restrained  device  which  is  suitable  for  Mossbauer  detection, 
that  is,  a  relative  velocity  output.  If  (65)  is  solved  for  a  step  in¬ 
put  in  acceleration,  A0,  the  following  solutions  in  x(t)  are  obtained: 


x(t)  =  A0 

t M) 1/2 

TT 

sin 

(  K  .1//2  ,  for  C=  0 

(68) 

x(t)  =  0, 

for  0  < 

C2 

<  i 

(69) 

ZKM 

x  • 

ft- 

II 

,  for 

C2 

<i 

(70) 

tKM 

The  condition,  C  =  0, 

in  (68)  is  difficult  to  achieve  in  practice. 

In 

order  to  achieve  high  sensitivity  in  this  system,  M  should  be  much 
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greater  thank  K;  however,  this  results  in  a  long  period  of  vibration, 

2  n  /  M  >  thus  limiting  the  time  response  of  the  system.  Final- 

ly  the  resolution  of  the  system  will,  in  all  probability,  be  defined 

by  uncertainty  in  K  and  not  by  the  pickoff  mechanism.  (69)  yields 

no  output.  (70)  presents  certain  practical  difficulties,  since  for  a 

large  A0)  x  may  rapidly  exceed  the  physical  dimensions  of  the 

instrument.  The  degree  of  resolution  that  may  be  attained  is  in- 

o 

herently  limited  by  the  condition  of  C  >4KM,  because  the  output 
is  proportional  to  M/C.  And  again,  performance  is  limited  by 
uncertainty  in  the  value  of  a  material  parameter,  in  this  case,  C. 

In  general,  simply  restrained  devices  are  not  widely  accepted 
as  the  best  approach  to  inertial  quality  accelerometers  (Ref  7:99, 
24:74).  Since  these  devices  are  limited  by  uncertainties  in  restraint 
parameters  such  as  C  and  K,  there  is  no  promise  that  simply  re¬ 
strained  accelerometers  can  be  improved  in  resolution  or  threshold 
sensitivity  by  the  application  of  a  sensitive  pickoff  system  employing 
the  Mossbauer  Effect. 

Application  to  Force  Balance  Systems 

Force  balance  accelerometers  offer  certain  advantages  over 
simply  restrained  devices  in  that  the  inertial  element  is  constrained 
to  small  displacements  about  a  null  position,  thus  minimizing  the 
non  linearity  and  hysterisis  inherent  in  end  range  inputs  to  an  elastic 
or  friction  constraint.  This  is  illustrated  by  the  Kearfott  Pendulous 
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Accelerometer,  Type  No.  F2401,  which  has  a  resolution  of  10^ 
over  a  range  from  5  x  10"^ g  to  20  g  (Ref  27:76-77),  and  by  the 
Litton  Model  A-200  Accelerometer  which  has  a  resolution  of  10^ 
over  a  range  from  10‘^  g  to  60  g  (Ref  27:92-93).  To  some  extent 
performance  of  a  force  balance  system  may  be  limited  by  such 
factors  as  electromagnetic  hysterisis  and  non  linearities  in  the 
forcing  device.  A  fair  indication  of  the  state  of  the  art  may  be 
obtained  by  considering  the  pendulous  integrating  gyro  accelerometer, 
abbreviated  PIGA.  The  PIGA  finds  wide  application  as  an  inertial 
quality  accelerometer  (Ref  7:99,  24:77). 

Pendulous  Integrating  Gyro  Accelerometer.  The  essential 
elements  of  the  PIGA  are  illustrated  in  Figure  13.  The  spinning 
wheel  of  a  gyro  creates  an  angular  momentum  vector  R,  given  by 

R=Iccspin  (71) 

where  I  is  the  moment  of  inertia  of  the  gyro  wheel.  The  vector  R 
is  constrained  to  move  about  two  axes,  indicated  by  an  inner  and 
outer  gimbal.  If  a  torque,  T  applied,  is  applied  to  the  inner  gimbal 
by  the  action  of  system  acceleration,  A  ,  on  the  two  masses,  M/2, 
then  the  system  will  precess  in  the  direction  of  uJ.  This  is  deduced 
from  the  law  of  conservation  of  angular  momentum,  expressed  in 
this  case  as 


where 


T  applied  =  cc  x  H 
T  applied  =  MAI  cos  (3 


(72) 

(73) 
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Figure  13 

Schematic  Diagram  of  a  Pendulous 


Integrating  Gyro  Accelerometer 
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p  is  the  angle  between  the  Spin  Reference  Axis  (SRA),  which  is  per¬ 
pendicular  to  co,  and  the  direction  of  TT  .  In  practice  p  is  kept 
very  small  (to  within  a  few  seconds  of  arc)  by  the  action  of  a 
servomechanism  which  turns  the  outer  gimbal  in  the  direction  of  co  . 

If  p  is  nearly  zero,  then  by  (72)  and  (73) 

co  =  MIA  (74) 

“TT" 

The  output  co,  is  directly  proportional  to  input  acceleration,  A,  over 
a  very  wide  range  of  performance.  In  practice,  the  output  is  usually 
taken  as  angular  position,  0  (  =  cot),  which  is  directly  proportional 
to  the  input  velocity.  Given  the  best  possible  design,  the  major 
uncertainty  due  to  material  properties  is  the  uncertainty  in  mass 
unbalance  (Ref  7:101).  In  the  case  of  an  inertial  quality  gyro  the 
uncertainty  would  be  about  15  x  10"®  cm,  expressed  as  uncertainty 
in  the  arm  length,  1  (Ref  7:103).  Even  for  an  arm  length  of  only 
one  centimeter,  this  respresents  a  small  inherent  error  of  only  about 

7 

one  part  in  10  . 

There  are  two  possibilities  for  using  the  Mossbauer  Effect 
in  conjunction  with  the  PIGA.  The  first  way  would  be  to  use  it  as 
a  position  sersor  for  detecting  p.  The  second  way  would  be  as  a 
means  of  measuring  cu,  the  output  angular  velocity. 

The  discussion  of  the  Mossbauer  Effect  as  a  position  sensor 
is  contained  in  Chapter  IV.  The  rotation  of  the  gyro  wheel  is  an 
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available  means  of  establishing  source  absorber  relative  velocity. 
co  spin  is  usually  high,  of  the  order  of  10^  radians/sec  (Ref  7:101, 
102).  For  an  absorber  mounted  just  one  centimeter  from  the 
center  of  rotation  of  the  gyro  wheel,  and  a  source  of  photons  aligned 
parallel  to  the  SRA,  by  (64)  an  angle  (3  of  .  02  seconds  of  arc  would 
result  in  a  relative  velocity  of  10“2  cm/sec.  However,  the  limita¬ 
tion  will  not  be  the  angular  speed  of  rotation  or  the  dimensions  of 
the  wheel  but  the  size  of  the  slit  in  the  shielding  that  must  be  employ¬ 
ed.  The  example  in  Chapter  IV  suggests  a  limitation  of  2"  of  arc. 

A.  pickoff  accuracy  of  2  seconds  of  arc  is  not  far  (if  at  all)  removed 
from  the  present  state  of  the  art.  For  instance,  the  M.I.T. 
Instrumentation  Laboratory,  Type  H  Integrating  Gyro  has  a  pickoff 
uncertainty  of  4  seconds  of  arc  (Ref  27:66,67). 

The  uncertainty  in  co  that  would  be  caused  by  a  mass  unbalance 
of  15  x  10“®  cm  is  of  the  order  of  10" ?  radians/second  (Ref  7:103). 

As  mentioned  before,  present  techniques  call  for  measuring  the  angu¬ 
lar  displacement,  o,  and  thus  determining  the  integral  of  accelera¬ 
tion,  which  is  velocity.  However,  any  uncertainty  in  co  will  be 
integrated  over  time  to  yield  an  uncertainty  in  0,  thus  an  error  in 
the  velocity  measurement.  It  is  instructive  to  estimate  whether, 
by  using  the  Mossbauer  Effect,  it  is  possible  to  measure  an  angular 
velocity  to  an  accuracy  commensurate  with  the  uncertainty  in  co  due 
to  mass  unbalance.  Typical  PIGA  values  are: 
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H  =  106  gm  cm^/sec 
M  =  500  gms 

(Ref  7:101-102).  Using  1=2  cm  in  (74),  u>=10-3a.  If  A  =  lg  = 

103  cm/sec2,  then  a>  =  l  rad/sec.  From  Figure  10,  10-6  rad/sec 
is  the  minimum  discernable  angular  rate.  This  corresponds  to  an 
uncertainty  in  acceleration  of  10 -3  cm/sec^  or  10"6  g,  and  a  resolu¬ 
tion  at  lg  of  106.  This  implies  measurement  of  relative  velocity  of 
10‘5  cm/sec  which  requires  a  minimum  counting  time  of  1  second. 

In  other  words,  the  system  will  have  a  comparatively  slow  response 
time.  By  the  conventional  method,  measuring  Q,  just  100  seconds 
of  system  operation  will  result  in  an  angular  displacement  error 
(AG  =  A co t)  of  10_5  radians  or  2  seconds  of  arc,  (If  A  o>  =  10“7 
radians/sec  due  to  mass  unbalance).  This  will  be  within  the  region 
of  uncertainty  of  the  pickoff  and  after  any  period  of  operation  con¬ 
siderably  longer  than  100  seconds  the  system  resolution  will  be 
limited  by  the  inherent  uncertainty  in  mass  unbalance  rather  than 
by  the  uncertainty  in  pickoff  accuracy.  The  conclusion  is  that  there 
is  no  advantage  to  be  gained  by  applying  the  Mossbauer  Effect  to 
measurement  of  the  output  angular  rate  of  a  PIGA. 

It  is  hoped  that  the  preceeding  discussion  has  revealed  to 
some  extent  the  orders  of  magnitude  and  degrees  of  tolerance  involved 
in  improving  the  performance  of  accelerometers  used  for  inertial 
guidance.  It  is  not  inferred  that  the  Mossbauer  Effect  cannot  be  used 
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to  improve  the  performance  of  any  type  of  force  balance  accelero¬ 
meter.  One  of  the  difficulties  in  seeking  to  apply  the  Mossbauer 
Effect  to  the  PIGA  is  the  inherently  low  sensitivity  of  w  to  A  in 
such  a  device.  This  may  be  inferred  from  the  relationship  u>=10-3a. 
In  this  case  it  is  necessary  to  accurately  measure  very  small  values 
of  a;  in  order  to  obtain  a  high  resolution  in  A.  One  can  conceive 
of  a  device  which  offers  a  more  favorable  relationship  of  co  to  A. 

Such  a  device  is  illustrated  schematically  in  Figure  14  and  will  be 
referred  to  as  a  centrifugal  reaction  accelerometer. 

Centrifugal  Reaction  Accelerometer.  Figure  14a  is  a  diagram 
of  the  lumped  mechanical  components  of  a  centrifugal  reaction  ac¬ 
celerometer.  The  accelerated  system  consists  of  an  axis  about 
which  is  suspended  a  mass,  M,  at  a  distance,  1,  by  a  weightless 
arm  which  makes  an  angle  (3  with  the  axis.  The  mass,  arm,  and 
axis  are  all  rotating  at  an  angular  rate  u>.  r  is  the  distance  from 
the  mass  perpendicular  to  the  axis.  The  mass  is  acted  on  by  two 
D'Alembert  forces,  due  to  the  system  acceleration  and  the  centri¬ 
petal  acceleration.  Taking  moments  about  the  pivot  joining  the  arm 
and  the  axis, 

M  r  uA  cos  p  =  MAI  sin  p  (75) 

Since  r  =  1  sin  p,  it  can  be  shown  that 

u?  =  A  / 1  cos  p  (76) 

Since  1  is  a  fixed  quantity,  when  u)  and  p  are  known,  A  is  determined. 
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Figui’e  14 

Schematic  Diagrams  of  a  Centrifugal  Reaction  Accelerometer 
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For  nominal  values  take  1  =  10  cm,  cos  p  =1.  Then  =  0. 1A. 

This  indicates  that  the  values  of  u>  in  the  centrifugal  reaction  ac¬ 
celerometer  will  be  correspondingly  larger  than  the  values  of  to 
in  the  PIGA. 

In  practice,  either  p  or  uj  may  be  held  to  a  null  value  through 
the  action  of  a  servomechanism,  or  both  p  and  to  may  be  measured. 
Figure  14b  shows  the  Mossbauer  Effect  used  to  measure  both  p  and 
co.  When  p  is  at  the  prescribed  value,  photons  from  the  source  pass 
parallel  to  the  axis  of  rotation  of  the  absorber  wheel.  Thus  there 
is  zero  relative  velocity  between  source  and  absorber  along  the  line 
of  photon  emission  and  a  maximum  of  Mossbauer  absorption  takes 
place.  When  p  is  not  at  the  proper  value,  the  slit  in  the  arm  will 
not  be  parallel  to  the  axis  of  rotation  of  the  absorber  wheel,  and 
there  will  be  relative  velocity  between  source  and  absorber  along  the 
line  of  emission.  may  be  measured  by  either  the  phase  difference 
method  or  the  amplitude  difference  method. 

The  resolution  that  may  be  obtained  is  estimated  in  the  manner 
previously  discussed.  The  size  of  the  slit  may  limit  resolution  in 
p  to  about  10-5  radians  or  2  seconds  of  arc.  If,  by  (76) 

A  =  1  cos  p  (77) 


then 

dA  =  -  1  sin  p  (78) 
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Dividing  (78)  by  (77),  and  making  the  approximation  A  A  =dA, 

A  A/  A  = -tan  p  A  p  (79) 

If  tanp  =  l  and  Ap  =10“^  radians,  then  A/  A  A  =  10^.  This  is  not 
exceptional  resolution  as  compared  to  the  PIGA.  It  should  be  noted 
that  by  making  tan  (3,  or  p,  very  small,  the  resolution  could  be 
increased.  However,  if  any  lateral  acceleration  is  present,  mak¬ 
ing  (3  small  will  lead  to  an  unstable  system.  This  is  inferred  by 
noting  that  as  p  becomes  very  small  the  moment  arm  for  lateral 
acceleration  becomes  large  compared  to  the  moment  arm  for  ac¬ 
celeration  along  the  axis  of  rotation. 

It  is  interesting  to  speculate  what  performance  could  be  achieved 
using  the  Mossbauer  Effect  to  measure  angular  velocity  if  unlimited 
resolution  in  p  were  achieved  using  some  means  other  than  the  Moss¬ 
bauer  Effect.  (77)  is  differentiated  as  follows: 

2  cc  dec  =  1  (80) 

d A  1  cos  p 

If  du)=Au>,  by  (80  and  (77) 

Aj  =  1  A  A  (81) 

~  ~T  ~~K~ 

From  Figure  10,  if  oc=  10  rad/sec,  A  /  AA  =  5  x  10^. 

This  is  promising,  however  it  must  be  accompanied  by  com¬ 
parable  resolution  in  p.  A  plot  of  resolution  in  acceleration  vs 
resolution  in  p  is  shown  in  Figure  15.  The  plot  assumes  tan  p  =1 
in  (79).  Unless  a  very  large  instrument  is  used,  no  known  method 
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will  provide  the  accuracy  of  angular  measurement  required  to  allow 
a  resolution  higher  than  10^  with  the  centrifugal  reaction  accelero¬ 
meter. 

One  further  point  is  the  question  of  whether  the  Mossbauer 
Effect  can  be  used  in  a  servo  loop.  When  it  is  considered  that  a 
relative  velocity  of  a  little  more  than  lO-^  cm/sec  will  destroy 
resonant  absorption  completely,  and  that  with  a  maximum  counting 
rate  of  10®/  sec  one  whole  millisecond  is  required  to  obtain  just 
100  counts,  it  is  difficult  to  justify  the  use  of  the  Mossbauer  Effect 
in  a  modern  control  system. 

If  the  examples  are  indicative,  then  the  application  of  the 
suggested  techniques  is  not  a  promising  method  for  improving  the 
resolution  of  force  balance  accelerometers.  However,  the  best 
that  can  be  said  is  that  the  discussion  emphasizes  the  fact  that  the 
resolution  of  such  accelerometers  is  limited  by  a  variety  of  factors. 
Even  assuming  a  technique  which  promises  extreme  accuracy  in  the 
measurement  of  a  particular  parameter  such  as  angular  velocity,  it 
is  often  useless,  considering  uncertainties  in  other  important  para¬ 
meters  such  as  mass  unbalance  and  angular  position,  to  apply  this 
technique  to  improve  the  accuracy  of  the  pickoff  system. 

In  the  light  of  these  considerations,  the  accelerometer  dis¬ 
cussed  in  the  following  section  is  of  particular  significance. 
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Application  to  Unrestrained  Devices 


An  interesting  example  of  an  unrestrained  accelerometer  is 
the  rotating  pendulum  accelerometer,  described  by  Samuel  Schal- 
kowsky  and  Henry  F.  Blazek  (Ref  25:469-473).  The  system  is 
illustrated  schematically  in  Figure  16.  A  mass  M  is  rigidly  sus¬ 
pended  from  a  pivot  0  at  a  distance  r.  The  mass  is  rotated 
through  an  angle  0  at  a  rate  to.  The  system  is  accelerated  at  the 
rate  A  in  a  direction  measured  from  the  y'  y  axis  by  the  angle  <)>. 
The  equation  of  motion  is 

I  d20  +  C  de  +  MAr  sin  (9  +  tj>)  =  Td  (82) 

"dt2  "dt 

where  I  is  the  moment  of  inertia  of  the  rotating  pendulum,  C  is 
the  coefficient  of  viscous  friction,  and  Td  is  the  driving  torque 
(Ref  25:471).  The  following  approximations  allow  an  analytic  solu¬ 
tion  of  the  equation: 

Let  0=o^ t  in  sin  (o+  cj>),  assuming  that  the  effect  of  the 
acceleration  produces  small  distortion  in  a  constant  angular  rate,  to  . 

Let  T^  =C  u&,  assuming  that  the  driving  torque  is  constant 
and  just  sufficient  to  overcome  the  effects  of  viscous  friction. 

The  approximate  solutions  of  the  equation  are: 

0  =coat  +  MrA  x  sin  ( u^t  +4>  +tan _1  C  ) 

I^a2  (  1+  C2  )  ,  /■;  (  I“a  ) 

(  T2^2) 

-  MrA  cos  ({>  (83) 

“a  c 
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Figure  16 

Schematic  Diagram  of  the  Rotating  Pendulum  Accelerometer 

(From  Ref  24:470) 


a 


Figure  17 

Phase  Difference  Method  Applied  to  the 
Rotating  Pendulum  Accelerometer  to  Measure 
Time  Increments 
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0  =  wa  +  _ Mr A. _ y_  X  cos  (  coat  +  4  +  tan'1  C  )  (84) 

Iwa  (  1  +  -C--T7-)  1/2  (  ^'a  ) 

(  &a2) 

(Ref  25:471).  It  is  now  noted  that  by  making  C  «  I  coa ,  the  solu¬ 
tion  in  0  becomes  independent  of  C;  that  is, 

b  =  aja  +  M  r  A  cos  (  u>a  t  +  cp )  (85) 

If  the  component  of  angular  rate  represented  by  the  last  term  can 
be  measured,  then  ooth  the  magnitude  and  direction  of  system  ac¬ 
celeration  can  be  determined:  Let 

0  =  coa  (86) 

where 

A  u ;a  =  Mr  A  cos  (  u>at  +  4. ) ;  (87) 

I  coa 

or,  considering  amplitude  only 

Atoa  =  Mr  A  (88) 

I  wa 

If  the  entire  rotating  system  consists  of  the  mass,  M,  suspended 

2 

at  the  distance  r,  then  I  =  M  r  ,  and 

A  o>a  =  A  (89) 

rwa 

It  is  desirable  to  make  ioa  as  small  as  possible  in  order  to  obtain 
a  high  sensitivity  of  Acoa  to  A.  If  o.'a  =  1  radian/ sec,  and  r=10  cm, 
and  if  A  cua  is  measured  to  10_D  radians  /  sec  using  the  Mossbauer 
Effect,  then  an  acceleration  of  10”^  cm/sec^  or  10"®  g  is  the  mini- 
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mum  detectable  acceleration,  or  threshold  sensitivity.  This  may 
be  reduced  even  further  by  reducing  wa.  The  limiting  condition 
is  that  C«Icoa. 

It  is  possible  to  apply  another  sensing  technique  to  the  rotat¬ 
ing  pendulum  accelerometer.  By  measuring  the  time  taken  by  the 
pendulum  to  pass  through  two  quadrants  and  comparing  it  to  the 
time  taken  by  the  pendulum  to  pass  through  the  remaining  two 
quadrants,  both  the  magnitude  and  direction  of  the  incident  accelera¬ 
tion  may  be  measured.  The  expressions  are 

A  =  Iwa3  (Ati2+Ato2)  (90) 

T MT 

4)  =  tan  -1  At2  /  Atj  (91) 

if  C  <<  Icoa  (Ref  25:472).  Atj  is  the  time  to  traverse  quadrants 
II  and  III  minus  the  time  to  traverse  quadrants  I  and  IV.  At2  is 
the  time  to  traverse  quadrants  I  and  II  minus  the  time  to  traverse 
quadrants  III  and  IV. 

Figure  17  represents  the  phase  difference  method  used  to 
instrument  the  rotating  pendulum  accelerometer.  Suppose  that  both 
absorbers  1  and  2  are  kept  in  phase  with  the  emitter  which  is  located 
adjacent  to  the  pendulous  mass.  Then  the  difference  in  phase  between 
the  excitation  voltages  supplying  the  coils  1  and  2  will  represent  a 
phase  change  due  to  Atj  or  At  2,  depending  on  the  orientation  of 
the  acceleration.  For  the  case,  At  2  =0,  a  maximum  phase  difference 
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would  be  equal  to  o>a  Atj.  Then  in  analogy  to  (34),  the  expression 
for  relative  velocity  between  source  and  absorber  becomes 

vr=Ve  sin  mat  -  Va  sin  cua  (t  +At^)  (92) 

This  can  be  approximated  (considering  absolute  values)  by 

vr  =Ve  coa  Atj  =Ru )a^  At  i  (93) 

Or,  equally 

Atl  =  (94) 

r^a 

By  using  (90)  and  (94),  an  estimate  can  be  made  of  the  minimum 

threshold  sensitivity  that  can  be  achieved  with  this  method.  The 

results  are  plotted  in  Figure  18,  which  also  shows  the  equivalent 

o 

results  from  (89).  In  both  methods  it  is  assumed  that  I  =  Mr  , 
vr  =  10-5  cm/  seC;  anc}  r  -  io  cm. 

Figure  18  indicates  that  extremely  low  values  of  threshold 
sensitivity  can  be  achieved  for  small  values  of  cua .  This  implies 
that  the  time  response  of  the  system  will  be  very  slow.  An  ap¬ 
plication  for  such  an  instrument  would  probably  be  to  instrument 
the  gravity  gradient  techniques  described  in  Appendix  A,  or  to 
monitor  the  performance  of  a  low  thrust  propulsion  device. 

The  rotating  pendulum  accelerometer  is  unique  in  that  by 
proper  design  the  output  can  be  made  independent  of  troublesome 
parameters  such  as  friction.  It  is  better  suited  for  an  application 
of  the  Mossbauer  Effect  than  most  accelerometers  because  the  slow 
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Threshold  Sensitivity  (g) 


Velocity  Measurement  -  eqn  (89) 

Phase  Measurement  -  eqns  (90)  and  (94) 


Figure  18 

Threshold  Sensitivity  of  a  Rotating  Pendulum  Accelerometer 
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period  of  rotation  should  afford  ample  time  for  counting.  In  a 
real  instrument,  of  course,  a  number  of  emitters  and  absorbers 
would  have  to  be  used  so  that  counting  time  would  be  maximized 
and  so  that  a  360°  profile  of  the  pendulum  velocity  could  be  obtained. 
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VI.  Conclusions  and  Recommendations 
The  purpose  of  this  paper  was  to  present  the  results  of  an 
investigation  to  determine  whether  the  Mossbauer  Effect  can  be 
applied  to  accelerometers  used  for  inertial  navigation.  The  investi¬ 
gation  has  not  led  to  a  final  conclusion,  but  it  is  believed  that  the 
material  presented  in  the  Thesis  supports  the  following  statements. 
Conclusions 

1.  The  Mossbauer  Effect  is  a  practical  means  of  elim¬ 
inating  recoil  energy  losses  from  the  emission  of  gamma  photons. 
When  these  recoil  losses  are  eliminated  the  photon  energy  can  be 
determined  to  an  uncertainty  limited  only  by  the  Heisenberg  Uncer¬ 
tainty  Principle.  This  allows  experiments  that  are  sensitive  to 

12  IS 

energy  differences  of  one  part  in  10  -10  ,  depending  on  the  pre¬ 

cision  of  the  experiment,  the  isotope  used,  and  the  time  available 
for  counting  photons. 

2.  Major  limitations  on  the  application  of  the  Mossbauer 

Effect  are  the  time  required  to  count  photons  and  the  energy  line- 

width  of  gamma  photons  from  available  isotopes.  The  use  of  the 
67 

isotope  Zn  ,  which  exhibits  the  smallest  linewidth,  is  impractical 
because  of  the  relatively  small  (1%)  fraction  of  zero  phonon  events 
that  occur  in  this  isotope  and  the  concomitant  small  cross  section 
for  resonant  absorption.  Fe®7,  which  has  a  Mossbauer  cross 
section  about  10^  that  of  Zn®7,  exhibits  the  second  smallest  energy 
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iinewidth  as  well  as  the  highest  fraction  of  zero  phonon  events  of 
any  Mossbauer  isotope.  These  factors,  together  with  the  desir¬ 
able  characteristic  of  a  Debye  Temperature  of  about  450° K  which 
permits  a  large  fraction  of  zero  phonon  events  at  high  temperatures, 
make  Fe^  the  most  suitable  isotope  for  a  practical  application. 

3.  Of  all  the  physical  phenomena  which  alter  Mossbauer 
absorption,  a  relative  velocity  between  source  and  absorber  pro- 
duces  the  most  pronounced  effect.  If  10°  photons  can  be  counted 
in  one  second,  then  a  relative  velocity  of  the  order  10"^  cm/sec 
can  be  detected  by  the  use  of  the  isotope  Fe^7.  An  acceleration 
of  3g  would  be  required  to  produce  the  same  effect  as  a  relative 
velocity  of  10“5  Cm/sec.  It  is  therefore  desirable  to  attempt  the 
use  of  the  relative  velocity  effect  rather  than  the  direct  measure¬ 
ment  of  acceleration  for  an  application  to  accelerometers. 

4.  There  are  several  possible  techniques  by  which  the 
Mossbauer  Effect  can  be  used  to  measure  the  velocity  of  an  instru¬ 
ment  component.  It  is  conceivable  that  angular  velocities  may  be 
measured  to  an  accuracy  of  10_6  radians/second  or  less.  There 
is  also  a  possibility  that  the  Mossbauer  Effect  can  be  used  to 
measure  angular  position,  but  accuracy  at  distances  of  10  cm  is 
limited  to  about  2"  of  arc  by  the  requirement  for  collimating  photons. 

5.  The  application  of  the  aforementioned  techniques  to 
accelerometers  is  limited  primarily  by  uncertainties  present  in 
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accelerometers  from  causes  other  than  the  resolution  of  the  pickoff 
system.  Table  II  presents  rough  estimates  of  the  resolution  obtain¬ 
able  from  vai'ious  systems  and  the  limiting  factors. 

Table  II 

Resolution  of  Various  Types  of  Accelerometers 

Type  Accelerometer  Maximum  Resolution  Limitation 

Simply  Restrained  10*^  -  10^  Friction-Hysterisis 

PIGA  10?  Mass  Unbalance 

Centrifugal  Reaction  105  Angular  Measurement 

It  appears  that  none  of  the  accelerometers  discussed  can  be  im¬ 
proved  in  resolution  by  application  of  the  suggested  techniques. 

6.  On  the  other  hand,  it  should  be  possible  to  measure 
extremely  small  values  of  acceleration  by  applying  the  Mossbauer 
Effect  to  the  measurement  of  angular  velocity  or  phase  difference 
from  a  rotating  pendulum  accelerometer.  The  threshold  sensitivity 
that  can  be  achieved  in  this  manner  is  of  the  order  10"®  g,  if  2  n 
seconds  is  the  period  of  the  rotating  pendulum.  For  larger  periods 
the  threshold  sensitivity  decreases  linearly. 

Recommendations 

Further  study  in  the  application  of  the  Mossbauer  Effect  to 
accelerometers  and  other  precision  devices  is  indicated.  It  is 
certain  that  this  paper  has  not  "uncovered"  all  the  possibilities. 

It  would  seem  probable  that  an  experienced  engineer  could  contrive 
a  useful  device  from  a  phenomenon  as  sensitive  to  small  energy 
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changes  as  the  Mossbauer  Effect. 

Experimental  verification  is  required  for  the  velocity  and 
position  measurement  techniques  described  in  Chapter  IV.  If 
proven  feasible,  the  techniques  should  be  inherently  useful.  The 
research  involved  in  this  study  has  not  revealed  any  other  techniques 
that  purport  to  measure  small  angular  velocities  to  the  degree  of 
accuracy  indicated  by  Chapter  IV. 

It  should  be  worthwhile  to  further  study  the  application  to 
the  rotating  pendulum  accelerometer.  If  the  small  threshold  sensi¬ 
tivity  indicated  by  this  study  can  be  reached  in  practice,  then  this 
will  provide  a  needed  means  of  monitoring  low  thrust  propulsion 
devices  and  a  means  of  realizing  the  navigation  techniques  suggested 
in  Appendix  A. 

Final  Statement 

The  Mossbauer  Effect  is  similar  to  many  other  discoveries 
of  contemporary  science  in  that  it  contributes  knowledge  and  new 
techniques  to  nuclear  and  solid  state  physics.  Just  as  other  modern 
discoveries  are  leading  to  profound  engineering  applications  such  as 
the  transistor,  the  LASER,  and  electric  propulsion,  it  is  probable 
that  the  Mossbauer  Effect  will  eventually  play  a  significant  role 
in  the  development  of  now  unknown  devices.  It  is  our  belief  that 
any  effort  spent  studying  the  Mossbauer  Effect  from  an  engineering 
standpoint  will  not  be  wasted. 
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\ppendix  A 

Navigation  by  Gravity  Gradient  Techniques 
This  Appendix  contains  an  example  of  how  specific  forces 
due  to  a  gravity  gradient  may  act  on  a  space  vehicle.  The  sug¬ 
gestion  is  made  that,  if  an  instrument  capable  of  sensing  these 
forces  becomes  available,  then  such  an  instrument  can  be  used 
effectively  as  a  means  of  navigation  and  control. 

In  Figure  19a  a  space  vehicle  is  represented  in  free  fall 
toward  the  center  of  attraction,  G,  of  the  earth.  At  the  center  of 
gravity  (C.G.)  of  the  vehicle,  the  kinematic  acceleration  of  falling 
just  balances  g.  However,  the  two  masses,  Mi  and  M2,  are 
subject  to  the  horizontal  component  of  g,  g  cos  p,  and  since  there 
is  no  horizontal  component  of  acceleration,  the  springs  must  supply 
a  tensile  force  to  restrain  the  masses.  The  magnitude  of  the  force 
is  small,  since  cos  p  is  L  /  r  where  L  is  the  distance  of  either  mass 
from  the  C.G. ,  and  r  is  the  distance  from  the  earth  center  to  the 
C.G.  A  nominal  value  of  the  specific  force,  for  L  =  1  foot,  r  = 

O 

earth  radius,  would  be  5  x  10"°  g. 

In  Figure  19b  another  case  of  specific  force  due  to  a  gravity 
gradient  is  represented.  Suppose  that  Mi  and  M2  are  oriented  along 
the  line  of  free  fall  as  shown.  Since  M*  and  M2  are  constrained 
by  the  structure  of  the  vehicle  they  must  fall  at  the  same  rate  as 
the  C.G.  However,  the  C.G.  accelerates  in  response  to  g,  whereas 
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the  specific  forces  acting  on  the  masses  are  gj  and  g2,  where, 
for  example 

gl  =  g  (  r  )  2  (95) 

(  r'  +  'L  ) 

in  accordance  with  the  well  known  inverse  square  law  of  gravity. 
Therefore,  the  spring  supporting  must  supply  a  compressive 
force  of  Mj  (g  -  gj).  The  nominal  value  of  such  a  force  is  small. 
For  L  =  1  foot,  r  -  earth  radius,  the  specific  force  would  be  about 

10"7  g. 

The  two  cases  just  discussed  can  be  generalized  for  an  ar¬ 
bitrary  orientation  of  the  vehicle.  In  Figure  19c  a  is  the  angle 
between  the  line  joining  Earth  center  and  the  C.G.  and  the  line 
Mj  and  M2.  A  good  approximation  for  the  specific  force  in  the 
springs,  Ag.  is 

Ag  =  -  (  2gs  r2  )  (  L  )  (1-3  cos2  a  )  (96) 

(  T2  )  (~T) 

where  gs  is  the  specific  force  of  gravity  at  the  Earth's  surface, 
and  R  is  the  radius  of  the  Earth.  A  tensile  force  is  taken  as  posi¬ 
tive  (Ref  2:46).  Using  (96)  one  can  determine  that  there  is  zero 
specific  force  for  four  vehicle  orientations: 

a  =  25°  16',  54°  44',  234°  44',  305°  16’. 

Since  specific  force  is  equivalent  to  acceleration  in  its  effect 
on  a  mass,  it  is  apparent  that  very  sensitive  accelerometers  could 
be  used  to  determine  the  orientation  of  the  vehicle  with  respect  to 
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the  attracting  body.  If  the  vehicle  orientation  is  compared  to  an 
inertial  reference  established  by  gyroscopes,  then  the  angular  rate 
of  the  space  vehicle  about  a  planet  can  be  determined.  The  mag¬ 
nitude  of  the  specific  forces  will  provide  information  about  the 
distance  of  the  vehicle  from  the  central  body.  This  is  summarized 
by  the  following  quotation  from  a  more  detailed  discussion  of  the 
same  topic  presented  here; 

"First,  it  would  be  possible  to  slave  an  in¬ 
ertial  platform  in  a  satellite  vehicle  into  continuous 
alignment  with  the  equipotential  surfaces  of  gravity 
and  thus,  by  purely  inertial  means,  to  determine 
continuously  the  attitude  of  the  vehicle  with  respect 
to  the  local  vertical.  Furthermore,  this  slaving 
could  also  be  applied  to  the  orbit  plane,  and, 
finally,  the  angular  velocity  of  the  satellite  in  its 
orbit  can  be  directly  derived. 

"It  follows  that,  if  the  gravitational  constant 
of  the  attracting  body  is  known,  the  radial  distance 
from  the  center  of  attraction  can  be  determined 
and,  conversely,  if  the  radial  distance  is  known, 
then  the  gravitational  constant  can  be  determined. " 

(Ref  2:45) 
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